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INTERPLANETARY MAGNETIC FIELD AND THE AURORAL ZONES* 


J. W. DungeyT 
Ionosphere Research Laboratory, The Pennsylvania State University, University Park, Pennsylvania 
(Received November 10, 1960; revised manuscript received December 22, 1960) 


The discovery’ of a regular interplanetary mag- 
netic field by Pioneer V has reawakened interest 
in Hoyle’s* suggestion that the primary auroral 
particles are accelerated at neutral points in the 
combination of an interplanetary field and the 
geomagnetic field. Hoyle pointed out that the 
latitude of the aurora would depend on the dis- 
tance of the neutral points from the earth and 
hence on the interplanetary field strength in the 
observed sense. The estimated particle energy 
was also reasonable. Dungey® discusses the ac- 
celerating mechanism. Here a qualitative model 
of the whole field is outlined and is found to be 
confirmed by the observed Sp current system. 

Consider a model with interplanetary plasma 
moving relative to the earth, this “wind” lying 
approximately in the ecliptic plane, and an inter- 
planetary field pointing roughly southward. The 
problem with no interplanetary magnetic field 
has been studied by several authors,‘ but the 


inclusion of this field alters the problem radically. 


The basic and awkward problem is that of the 
flow of plasma round the earth. This has not 
been solved but will be sketched using the physi- 
cal picture of hydromagnetics. The flow ina 
plane containing the neutral points is shown in 
Fig. 1. The flow near the neutral points is con- 
trolled by the strong current density existing 
there (3). The reverse flow between the neutral 
points is important; it has to flow round the 
earth, but does so in a normal aerodynamic way. 
A steady laminar flow will be assumed here for 
simplicity, but it should be noted that large var- 
iations of the field were detected by Pioneer I.° 


Denoting the local wind velocity by u, the electric 
field E is approximately -UxH/c everywhere 
outside the ionosphere except near the neutral 
points. In a steady state E has a potential, which 
is constant on a line of force, with the same ex- 
ceptions. Far from the earth E points out of the 
paper and the direction is still much the same 
near the neutral points and in the equatorial part 
of the region of reverse flow. In order to deduce 
the ionospheric currents, the topology of the 
magnetic field must be considered. 

In Fig. 1 there will be two lines (not in the 
plane shown) connecting the neutral points and 
together forming an approximately circular 
closed curve C near the equatorial plane. The 
lines from one neutral point will cover a surface, 
topologically similar to half a cylinder, extend- 





-— LINE OF FORCE 
==> DIRECTION OF FLOW 


FIG. 1. Interplanetary plasma flow in a plane con- 
taining neutral points. 
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ing from C to infinity to the south and a surface 
like half a doughnut bounded by C and intersect- 
ing the earth’s surface in a closed curve A in the 
northern hemisphere. The lines from the other 
neutral point complete the cylinder and doughnut 
and each neutral point has two lines on the other’s 
surfaces. The cylinder and doughnut separate 
regions where lines of force have both, one, or 
no feet on the ground, the other ends going to 
infinity. 

The plane of Fig. 1 is at constant potential by 
symmetry. Consider lines of force leaving the 
earth in a circle of longitude perpendicular to 
Fig. 1. The electric field fed in near the pole 
points out of the paper. Lines outside the curve 
A (that is, at lower latitude), however, pass 
through the equatorial plane and E here is out of 
the paper; following the potential up lines of 
force then gives a reversal of E at the curve A , 
because the lines from the equatorial plane have 
turned through more than 90°. Consequently the 
equipotentials in the ionosphere are as shown in 
Fig. 2. These equipotentials resemble the Sp 
current lines and the explanation is that Hall 
conductivity dominates in the ionosphere. Study 
of the conductivities o,, 0,, and o, in the iono- 
sphere shows that the electric field will extend 
with only slight diminution down to the E layer, 
so that the integrated conductivity is appropriate 
and the current density is nearly perpendicular 
to E. The phase of the current may be obtained 
by taking the direction of jx at the pole to be 
opposite to that of ixH far out and then a wind 





Wind direction } 


FIG. 2. Equipotentials in northern hemisphere for 
plasma winds of Fig. 1. 


from the sun fits the observations. 

The fit with observed Sp suggests the truth of 
the model. Sp must involve currents in the iono- 
sphere driven from outside and the sudden rever- 
sal at the auroral zone seems to require neutral 
points. If the interplanetary field is northward 
the topology is quite different, the surface of 
lines of force enclosing the earth but not meeting 
it. It is therefore concluded that there is an ap- 
proximately southward field. 

One consequence of the flow in Fig. 1 is a mix- 
ing of the interplanetary plasma into the outer 
atmosphere, which will upset the equilibrium 
under gravity which would exist in the absence 
of such a flow. The whistler data suggest that 
plasma is being lost far out.® 

The connection between the neutral points and 
auroras is obvious in this model, but it remains 
to study the motion of individual particles and 
the effect of turbulence. The outer Van Allen 
belt also remains to be discussed; the return 
flow in Fig. 1 must not upset its symmetry about 
the geomagnetic axis. The model predicts an 
asymmetry for auroras: for a wind from the sun 
there should be proton auroras before midnight 
and electron auroras after midnight. Recent re- 
sults’ from IGY data are of great interest in 
connection with the model presented here and 
promise substantial advances in our understand- 
ing of the subject. 





*The research reported in this work has been spon- 
sored by the Geophysics Research Directorate of the 
Air Force Cambridge Research Center, Air Research 
and Development Command, and, in part, by the Na- 
tional Science Foundation. 

TPresent address: Atomic Weapons Research Estab- 
lishment, Aldermaston, England. 
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COSMIC -RAY NEUTRON INCREASE FROM A FLARE ON THE FAR SIDE OF THE SUN 


H. Carmichael and J. F. Steljes 
Atomic Energy of Canada Limited, Chalk River, Canada 


D. C. Rose 
Division of Pure Physics, National Research Council, Ottawa, Canada 


and 


B. G. Wilson 
Sulphur Mountain Laboratory, University of Alberta, Calgary, Canada 
(Received December 27, 1960) 


Two conspicuous bursts of cosmic-ray neutrons 
produced by large visible solar flares were ob- 
served at many laboratories throughout the world 
on November 12 and 15, 1960, and these events 
are now being intensively studied. The purpose 
of this note is to draw attention to a third and 
rather inconspicuous injection of solar cosmic 
radiation that occurred on November 20, 1960, 
and to point out the interesting fact that the radia- 
tion in this case seems to have originated in a 
large unseen flare which occurred at a solar lon- 
gitude some 120°W of the center of the solar disk; 
in fact, about 30° behind the west limb of the sun. 
The November 20 injection was noted in Canada 
at Deep River and at Sulphur Mountain and has 
been found to be present also in the Ottawa and 
Churchill neutron monitor data. Results from 
Resolute are not yet available. 

The increase started at 2055+ 10 U.T. on No- 
vember 20. The preliminary report of solar 
activity’ does not list any large flare such as in 
the past has always been seen at the time of a sea 
level solar cosmic-ray increase, while the list- 
ing of numerous small flares shows that observ- 
ing conditions were good. The most likely place 
for a large flare to have occurred was McMath 
plage region 5925 which had already produced 
major flares (of importance 3 and 3+) on No- 
vember 10, 12, and 15. This region was still 
active as it passed round the west limb with 
loops, surges, and a succession of small flares. 
Assuming a 27-day rotation period and taking the 
3+ flare at 1322 U.T. November 12 at 10°W of the 
center of the solar disk,’ the corresponding posi- 
tion for 2055 U.T. on November 20 is 120°W and 
we postulate that a major solar flare occurred 
there about that time. 

This postulate is strongly supported by the ob- 
servations of solar radio noise made by Coving- 
ton® and analyzed in the accompanying Letter.* 

Detection of solar radio noise requires line-of- 
sight propagation from the source to the observ- 


ing station and consequently noise outbursts are 
observed with comparable frequency from both 
the east and the west half of the solar disk. In 
the case of solar cosmic rays, however, it has 
been noted that solar cosmic radiation seldom 
reaches the earth from flares near the center 
and on the eastern half of the disk. If it does, 

the signal shows characteristically a very slow 
rise. From flares on the west half the increases, 
whenever they occur (only 10 solar increases in 
all at sea level prior to the one being discussed 
have been recognized in the past 20 years), are 
fast rising and tend to start sooner in relation to 
the time of the visible flare. This behavior has 
recently been® attributed to the existence of con- 
figurations of magnetic lines of force above a 
strongly active solar region which stretch in an 
easterly curving path, as a result of the rotation 
of the sun, from the sun to beyond the orbit of the 
earth. In the present case we must assume that 
these magnetic lines of force formed an easy path 
to the earth for high-energy protons generated 
some way behind the western limb of the sun. 

The evidence for a neutron increase due to 
solar cosmic rays is exhibited in Fig. 1. The 
pressure-corrected hourly totals measured by 
the standard neutron monitor at Deep River cover 
a period of 15 days immediately following the 
solar cosmic-ray increase of November 15. Dur- 
ing this period the intensity is recovering from 
the Forbush decreases that occurred on Novem- 
ber 12, 13, and 15. Large diurnal and other slow 
fluctuations are apparent. There is also a small 
Forbush decrease at 0000 U.T. on November 25. 
At 2100 hours Universal Time on November 20 
(Fig. 1) a sudden increase of neutrons is to be 
seen at Sulphur Mountain, Deep River (Ottawa is 
similar), and at Churchill, but the mesons at 
Sulphur Mountain did not exhibit a corresponding® 
increase. It is also evident from Fig. 1 that the 
sudden neutron increase at Sulphur Mountain 
(altitude 2283 m), relative to the normal count 
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FIG. 1. Pressure-corrected hourly totals from 
standard neutron monitors and from a meson telescope. 
The ordinate scale divisions are at intervals of approxi- 
mately 5% for the neutron monitor graphs and 2.4% 
for the meson telescope graph. 


rate there, was more than double the size of the 
increase seen at the sea level stations. The char- 
acteristic shape of this increase, the fact that it 
was relatively much larger at a mountain station, 
and its absence in a meson telescope at that same 
station, indicate that it was produced by solar 
protons. 

The Sulphur Mountain and Deep River neutron 
totals are plotted on an expanded time scale in 
Fig. 2 to show the shape of the increase. For 
this figure the neutron counting rate at Deep River 
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FIG. 2. Uncorrected 10-minute totals for Sulphur 
Mountain and pressure-corrected 5-minute totals for 
Deep River. The basic (i.e., 100%) rates for these 
graphs are 100600 per hour at Sulphur Mountain and 
149 000 per hour at Deep River. 


has been more than doubled by including the counts 
from a graphite-moderated neutron monitor. 

From Fig. 2 the onset of the increase appears to 
have been at 2055+ 10 U.T. and the rise to maxi- 
mum took about one hour. 





‘Billings, Trotter, and Wetmore, High Altitude Ob- 
servatory Preliminary Report on Solar Activity, 
Boulder, Colorado, 1960 (unpublished). 

"Helen Dodson-Prince, McMath-Hulbert Observatory 
(private communication, November, 1960). 

3Daily values of solar flux at 2800 Mc/sec (10.7 cm) 
and outstanding events recorded at National Research 
Council, Ottawa, Canada, in November 1960; Decem- 
ber, 1960. 

‘a, E. Covington and G. A. Harvey, following Letter 
[Phys. Rev. Letters 6, 51 (1960). 

5k, G. McCracken and R. A. R. Palmeira, J. 
Geophys. Research 65, 2673 (1960); and references 
therein. 

Note that the vertical percentage scale on which the 
meson totals are plotted in Fig. 1 has been expanded 
with respect to that used for the neutron graphs by a 
factor 2.1 to make the general variations look the same. 
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10.7-cm SOLAR NOISE BURST OF NOVEMBER 20, 1960 


Arthur E. Covington and Gladys A. Harvey 
Radio and Electrical Engineering Division, National Research Council, Ottawa, Canada 
(Received December 27, 1960) 


The occurrence of a flare upon the far side of 
the sun as described in the previous Letter by 
Carmichael et al.’ is an extremely likely event 
as inferred from analysis of two sets of 10.7-cm 
solar noise observations taken at the Radio Ob- 
servatory of the National Research Council, 
Ottawa, Canada. 

The first set of observations was made with a 
1.3-meter parabolic reflector and a radiometer 
which records the solar emission from the whole 
disk from sunrise to sunset. A tracing of the per- 
tinent portion of the record for November 20, 
1960 from 19:00 U.T. to 22:00 U.T., is shown in 
Fig. 1. The burst is the enhancement above the 
daily quiet level—shown by the dotted line and 
estimated from the level on the previous day. 
The burst features from 20:30 until sunset at 
21:20 have been modified by the interference 
pattern formed by the direct solar ray and the 
ray reflected from the earth. The times of start 
of the main burst features are given in Table I. 
The preliminary listing of flares* for this period 
has been examined and the times of start of the 
six flares reported in this period are also listed 
in Table I. The flares occurred in two regions— 
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FIG. 1. Record of 10.7-cm solar noise burst, 


November 20, 1960. 


one, the very active region of interest (HH 25); 
the other, a smaller region in the same latitude 
as HH 25 but following it by a few days (FF 32). 
Comparison of the times of start of the bursts 
and flares within the time interval customarily 
used for associating 10-cm bursts and flares 
show that two of the flares, 20:09 and 20:19, 








Table I. Comparison of flare and burst features. 
Flare features Burst features 
Beginning Importance Region Beginning Flux Remarks 
19:27 i- FF 32 No evidence of flux increase. 
19:39 7 Small rise and fall of flux. 
19:45 7 Small rise in flux to a steady level. 
19:55 1- FF 32 No evidence of flux increase. 
20:09 1- FF 32 20:09 11 Small rise and fall of flux on steady level. 
20:16 1- FF 32 No evidence of flux increase. 
20:17 1- HH 25 No evidence of flux increase. 
20:19 1 HH 25 20:19 7 Small rise in flux; no appreciable fall. 
20:23 400% Large increase in flux; continuing until 


end of observations. 





* Peak value obtained from low gain —high speed channel. Not illustrated. 
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coincide well with burst features. The lesser of 
these, at 20:09, occurs in the lesser and follow- 
ing region (FF 32) which is just inside the limb 
around which the main region (HH 25) is disap- 
pearing. The larger flare, at 20:19, occurs in 
the main region (HH 25) which at this time is 
mostly around the limb. The other small flares 
are not associated in time with any burst fea- 
tures— which is not unusual, in that this is the 
case for the majority of the smallest flares. The 
unusual feature of the burst at 20:23 is that this 
is a radio event of more than average intensity 
with a flux of 400x107? (watt/m*)/cps and a 
duration of over 45-minutes, while the only flare 
listing for approximately this time is of impor- 
tance 1 and duration 4 minutes. Comparison of 
10-cm events and flares over a long period of 
time has indicated that bursts of this magnitude 
are almost always associated with flares of im- 
portance 2 at least. This supports the suggestion 
of an intense flare on the far side of the sun oc- 
curring in the active region which has largely 
disappeared around the limb. Possibly the short- 
duration flare reported is a small fragment on 
the limb of the much larger hidden flare. 

A similar situation was observed previously on 
June 9, 1959. Analysis* showed in this case 
that the hidden flare occurred in a region which 
was just about to appear on the eastern limb. In 
this critical position, a direct line of sight is 
possible only for the radio emissions which are 
produced in elevated regions above the flare. It 
has also been suggested* that this event produced 
x-ray emissions. 

The second set of solar noise observations was 
made with a 185-meter aperture interferometer 
system with a fan-shaped beam 1’ E.W. by 2° 
N.S. The daily strip scans of the solar surface 
are made at 16:30 U.T. and enable the bright 
radio emissive regions to be accurately located 
in one dimension. On November 4, a radio 
emissive region appeared on the eastern limb, 
crossed the central meridian on November 12, 
and was last seen on the western limb on Novem- 
ber 20. The relation of the trajectory of the 
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radio emissive region with respect to the solar 
disk indicates that the region was 20° to 30° 
north of the solar equator. This corresponds to 
the location of the optical region HH 25. In addi- 
tion, the scan taken during the post-burst in- 
crease associated with the November 12 erup- 
tion clearly indicates increased emission from 
this region. Examination of the provisional sun- 
spot positions as reported by the U. S. Naval 
Observatory showed that the sunspots of this 
region were first seen near the eastern limb on 
November 6 and last seen on the western limb 
on November 18. The longer visibility of the 
radio region in comparison with the optical 
manifestations is typical® and shows again that 
the 10-cm radio emission is observed even when 
the optical center is located on the far side of 
the limb. This emission is from a hot coronal 
condensation located some 3% of the solar ra- 
dius above the photosphere. The surface area 
of the radio region seems related to bright 
plages with observed E. W. extents ranging from 
1 to 6 minutes of arc. 

When the reported sunspot positions are pro- 
jected to give a value for November 20, it is 
found that the region extends from 20° to 30° 
around the far side of the western limb. Since 
radio emission from the coronal condensation in 
this position was observed, it is reasonable to 
assume, in the absence of an intense flare on the 
visible disk, that a hidden intense flare was 
responsible for the burst of noise observed at 
20:23 U.T. 





'H, Carmichael, J. F. Steljes, D. C. Rose, and 
B. G. Wilson, preceding Letter [Phys. Rev. Letters 
6, 49 (1961)). 

Billings, Trotter, and Wetmore, High Altitude Ob- 
servatory Preliminary Report on Solar Activity, 
Boulder, Colorado, 1960 (unpublished). 

3R. J. Coates, A. E. Covington, and S. Edelson, 
Astron. J. 64, 326 (1959). 

‘J. Kleczek and L. Krivsky, Nature 186, 1035 (1960). 

5A. E. Covington and G. A. Harvey, Astrophys. J. 
132, 435 (1960). 
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HIGH-FREQUENCY POTENTIAL WELL FOR DENSE PLASMAS 


Albert J. Hatch 
Argonne National Laboratory, Argonne, Illinois 
(Received December 27, 1960) 


The stable dynamic containment of charged 
particles and low-density plasmas in nonhomo- 
geneous high-frequency electric or electromag- 
netic fields is now firmly established theoreti- 
cally'~* and experimentally.*~® In electromag- 
netic fields the containing force on electrons is 
derivable from a potential function® 


’ =(e?/4mw*)(E*) |. (1) 


Here e and m are the charge and mass of the 
electron, w is the radian frequency of the field, 
and 3(E*),, is the time-averaged square of the 
electric field (E- -E) ayy where E = =E exp(jwt). This 
scalar point function generates three-dimensional 
potential wells in the vicinity of the central field 
nodes in certain resonant-cavity standing wave 
modes such as the electric quadrupole mode. Be- 
cause of the necessity that the field penetrate the 
plasma in order to act on the individual particles, 
the plasma density is limited to values which 
make w,>w, where Wp =(ne*/me o)” is the char- 
acteristic radian frequency of the plasma and n 

is the electron density. Consequently, there has 
been only mild interest in the basic and applied 
aspects of @ potential wells. 

The purpose of this Letter is to describe a 
newly recognized absolute potential well which 
theoretically appears to be capable of contain- 
ing dense plasmas (wp> w), and thus does not 
suffer from the density limitation imposed on the 
@ well. The present analysis arose out of theo- 
retical considerations of the results of an exper- 
imental analog study of the equilibrium properties 
of simulated plasma cores in resonant cavity 
fields.’ 

If a small volume of sharply bounded dense 
neutral plasma is placed in a resonant cavity 
field, the penetration of the high-frequency £ and 
H fields into the highly conducting volume of 
plasma will be limited to a very small skin depth. 
It is therefore convenient to consider the interior 
volume of the plasma as field free and to assume 
that the plasma can be represented by a small 
rigid neutral conducting sphere. This sphere 
can be considered as a test volume whose suc- 
cessful containment by a particular cavity mode 
indicates the strong possibility of similar con- 


tainment by that same mode of a somewhat larger, 
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less well defined, nonrigid, nonspherical volume 
of dense plasma. 

The potential of a perfectly conducting sphere in 
a quasi-static parallel E field® is Ug =-77,*€,(E”) ay. 
Similarly, for a perfectly diamagnetic sphere in 
a quasi-static parallel H field the potential is 
Um = 377° HoH") ay. The net potential of this 
sphere in superposed parallel E and dH fields is 
therefore 


T= 07 (4H (H*), - €9(B*),.)- (2) 


This expression appears in the theory for the 
measurement of cavity fields by frequency per- 
turbation due to the presence of a small metallic 
sphere’ and therefore is not new in itself. The 
new step reported here involves the use of T as 

a field property with which the containment cap- 
abilities of various resonant -cavity field config- 
urations can be assessed. In order to apply 

Eq. (2) to the nonuniform fields in cavities it is 
necessary that 7,/A«1, where A=2mc/w, to in- 
sure that the local fields in which the sphere is 
placed be approximately parallel. The net aver- 
age body force acting on the sphere is -AT which 
phenomenologically is due to the dielectrophoretic 
and diamagnetic behavior of the sphere in the non- 
homogeneous field. Various aspects of the signif- 
icance of this force to the containment problem 
have been discussed recently by Yankov.’°® 

The exact calculation of T is greatly simplified 
by letting , approach 0, in which case T becomes 
a scalar point function of a point sphere. YT has 
been calculated for several spherical cavity modes 
for which the field components E(r, 6, ¢) and 
H(r, 6, ¢) are given in standard references," 
namely, electric dipole (TM,,,), magnetic dipole 
(TE,,.), electric quadrupole (TM,,,), magnetic 
quadrupole (7E,,,), and magnetic sextupole (TE,,,). 
Of these modes only the last two possess suitable 
absolute potential wells, and only the TE,,,, mode 
well is here discussed in detail. (A comprehen- 
sive report on the various modes is in prepara- 
tion.) 

The field structure of the TE,,, magnetic quad- 
rupole mode in the spherical cavity is shown 
qualitatively in Fig. 1. This is a cross section of 
the field in the (vy, 8) plane. The broken lines 
represent the H field, and the dots and crosses 
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FIG. 1. Field configuration of the TE.) magnetic 
quadrupole mode in a spherical cavity. Magnetic field 
lines are shown as broken lines; the dots and crosses 
represent the electric field. 


represent the E field whose lines are concentric 
circles centered on the polar axis. There isa 
node in both the E and H fields at the center of 
the cavity. 

Values of T(r, 6) for the magnetic quadrupole 
mode in a sphere of radius R are shown in Fig. 2 
for the upper right quadrant only. The radial 
variation of T has been plotted on an arbitrary 
scale at 10° increments of 8. In each curveT is 
zero at the center of the cavity where there is a 
field node, becomes positive for small values of 
r/R, and reaches a crest in the vicinity of 
0.20<7/R <0.43 as shown in Fig. 2. The min- 
imum crest height occurs at r/R =0.215, @=65°; 
the maximum crest height is on the polar axis, 
@=0°, at r/R =0.43. The crest at the equator, 
@=90°, also occurs at r/R =0.43. The ratio of 
the maximum to the minimum crest height is ap- 
proximately 9:1. The well thus formed inside 
the T crest is symmetrical in the four quadrants. 
The volume of the well inside the crest is 2.6% 
of the volume of the cavity; the volume of that 
portion of the well enclosed within the equipoten- 
tial drawn through the minimum crest point, as 
shown in Fig. 2, is 0.6% of the volume of the 
cavity. Neither of the configurations suggested 
by the crest or by the equipotential contours are 
necessarily indicative of the configuration of a 
finite dense plasma contained in the well. Indeed, 
the work of Berkowitz et al.” indicates that a 
more likely configuration is that of a hypocycloid 
of revolution. The apparent well in the vicinity of 
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FIG. 2. Potential function T(r, @) in the TE») mag- 
netic quadrupole mode in a spherical cavity. The po- 
tential well lies between the origin and the crest, and 
is symmetrical in the four quadrants. 


r/R =0.6 and 6=50° occurs at an E-field maximum 
rather than at a field node, hence is unsuitable 
for containment of dense plasmas. 

The necessary conditions for the existence of 
an absolute T well are that there be a point in the 
electromagnetic field where (a) T=E =H =0, 

(b) aT /ar=0, and (c) 8?T/a7*>0. Consideration 

of the field structure of the various modes in 
simple cavities (i.e., spherical, cylindrical, or 
rectangular) indicates that the above conditions 
exist simultaneously only in the magnetic quad- 
rupole and higher order magnetic multipole modes. 
It is significant that the magnetic quadrupole 

mode is the resonant cavity version of the static 
magnetic “cusp” field which has been shown” to 
have stable equilibrium properties for the con- 
tainment of dense plasmas. The present result i 
is also consistent with the levitation experiments 
in audio-frequency magnetic quadrupole fields 
described by Okress et al.,’* and provides an 
alternate explanation of the phenomena which they 
observed. 

The existence of the T wells (and hills) indi- 
cates the possibility of new basic phenomena in- 
volving the interaction of plasmas with electro- 
magnetic fields and suggests new approaches to 
the applied problem of the containment of dense 
plasmas. In either case one must go beyond the 
present analysis and consider additional factors 
such as the effects of coexisting @ potential wells 
or hills, high rf power requirements (of the order 
of kilowatts to produce pressures of dynes/cm’), 
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competition from other effects, and alteration 
of the T well by the finite size of a contained 
plasma core. In regard to this last point, con- 
sideration of the analog experiment’ (unpublished 
results) and the levitation experiment** indicates 
that an T well exists for test-sphere diameters 
up to about y,/R < 0.1 in the present case. 

It is a pleasure to thank A. Strecok for pro- 
gramming the computations, and G. R. Ringo 
and F. E. Throw for helpful comments. 





*Work performed under the auspices of the U. S. 
Atomic Energy Commission. 
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ASPHERICAL 3d ELECTRON DISTRIBUTION IN Nit+T 


Harvey A. Alperin 
U. S. Naval Ordnance Laboratory, Silver Spring, Maryland and 
Brookhaven National Laboratory, Upton, New York 
(Received December 22, 1960) 


The purpose of this Letter is to report the re- 
sults of the first measurement of an aspherical 
magnetic form factor of an antiferromagnetic 
single crystal, nickel oxide, using unpolarized 
neutrons. The simple structure, if this com- 
pound is purely ionic, makes it possible to ex- 
plain the resultant e. symmetry that is obtained. 
This differs from previous (polarized beam) 
measurements’ ~* where no simple explanation of 
the asymmetries is possible because the mate- 
rials studied were either metals or (in one case) 
had the complicated spinel structure. The form 
factor of Nit*+ is found to lie above the free- 
atom curve in contrast to the case of Mn**+ where 
the converse is true. These results are of inter- 
est because of the very recent Hartree-Fock spin 
polarized calculations for Ni** by Watson and 
Freeman.° 

Unpolarized neutrons were used since the more 
accurate polarized neutron beam technique® is 
not applicable to antiferromagnets like NiO that 
do not have superimposed nuclear and magnetic 
reflections. The single crystal used for this in- 


vestigation was cut in the shape of a cylinder 

(8 mm high, 2.3 mm in diameter) from a large 
boule grown by flame fusion. Integrated inten- 
sities of all nuclear and magnetic reflections 
were measured at room temperature up to a 
value of siné/) = 0.78 using a wavelength of 
1.046 A. From the heavily extinguished nuclear 
reflections a mosaic spread parameter of 112 
seconds of arc was determined. This resulted 
in secondary extinction corrections to the mag- 
netic intensities of only 12% for the strongest 
reflection (111), of only 3% for the next strong- 
est (311), and of a completely negligible amount 
for all the other weaker reflections which occur 
at higher Bragg angles. 

Effects due to double Bragg scattering became 
particularly noticeable when measuring the very 
weak outer reflections (signal to noise ratios in 
the range 0.05 to 0.1). In order to help correct 
for these effects, additional room temperature 
measurements were taken at a second wavelength 
of 1.005 A. By taking account of the changes that 
occurred in the peak shapes upon changing the 
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wavelength, it was possible to obtain estimates 
of the true intensity of the reflections under 
study. However, in many cases the error en- 
tailed in this procedure was rather large. 

The experimental values for the magnetic form 
factor are plotted in Fig. 1 as the circles. The 
values were put on an absolute scale by using 
the previous result’ that the value of the form 
factor for the (111) reflection is 0.94. The im- 
portant observation to be made is that values of 
the form factor are obtained which depart radi- 
cally from a smooth curve. In particular, outer 
reflections which occur at the same Bragg angles 
show very different values for their form factors. 
The existence of such a multivalued form factor 
demonstrates that the unpaired 3d electrons are 
distributed aspherically about each Ni** ion. 

To investigate the possibility (however unlikely) 
that the observed results might be due to the ef- 
fect of anisotropic motions, either of the atoms 
with respect to one another or of the atoms with 
respect to their electron clouds, the (11, 1, 1) 
and (5, 7, 7) reflections, which occur at the 
same Bragg angle, were measured after cooling 
the crystal to 117°K. If the effect were due to 
anisotropic vibrations one would expect the ratio 
of the form factors of the two reflections to ap- 
proach approximately unity as the temperature 
approaches 0°K. Within the limits of experi- 
mental error, this ratio did not change upon 
cooling from 297°K to 117°K. 


The aspherical unpaired 3d electron distribu- 
tion observed experimentally is just what is ex- 
pected due to the splitting of the energy levels of 
the d electrons in the field of the anions. Accord- 
ing to crystal field theory,® in an electric field 
of octahedral symmetry the five-fold degenerate 
levels split into a low-lying triplet of to, sym- 
metry and a higher doublet of €g Symmetry. 
Since the levels fill up according to Hund’s rule, 
the two unpaired electrons which give rise to the 
magnetic moment should then exhibit pure eg 
symmetry.® 

The form factor can be written as the sum of 
a spherical part and an aspherical part, f=/f5+/f,- 
The part f4 (having pure €g Symmetry in our 
case) can be calculated following the methods of 
Weiss and Freeman.*° However, in order to 
compare the experimental results directly with 
those predicted for an f, of eg symmetry, it is 
necessary first to determine fg. Since at the 
present time there is no adequate way of calcu- 
lating f; for atoms bound in a solid, fs was 
chosen to yield the best agreement between the 
experimental and calculated values of the form 
factor. It was found to be possible (also most 
convenient and instructive) to do this by scaling 
the free-atom spherical form factor" for Ni**, 
i.e., fg(x)=folkx), where fp is the free-atom form 
factor, x = sin@/r, and the scale factor k = 0.83. 
In Fig. 1, fg is plotted as the solid curve.’ The 
final calculated values of f are denoted by the 
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triangles. Reasonable agreement is thus obtained 
between the experimental form factor values and 
those computed assuming the unpaired electrons 
to be in a state of pure e, symmetry. 

It is important to observe that the scale factor, 
k, here determined experimentally, is such as 
to make the unpaired 3d charge distribution for 
Ni** much more compact in the solid than it is 
for the free atom. This is in contrast to the 
case of Mn** where experiments” show that the 
charge distribution is expanded in the solid. 

These experimental results may be compared 
with the recent calculations of Watson and Free- 
man’ for the Ni** ion. These Hartree-Fock self- 
consistent field calculations allow the wave func- 
tions of electrons with opposite spins to have 
different radial dependencies (spin polarization) 
and lead to a contraction of the charge distribu- 
tion (both for the free-atom case and the case 
where the Ni** ion is placed in an octahedral 
array of point charges).* Unfortunately, the 
magnitude of the contraction is much too small 
to explain the observations reported here. Never- 
theless, the important fact that the relation of 
the observed f; to the free-atom fg is just oppo- 
site for the cases of Nit*+ and Mn** would lead 
one to look for the origin of this effect in the out- 
standing difference between the two ions: namely 
their differing spin configurations. These ex- 
perimental results also serve to suggest that 
whereas effects such as spin polarization and 
crystalline environment have a large influence 
on fs, their effect on f4 is small. 

The author wishes to thank Professor E. Uchida 
and Dr. Y. Nakazumi for the single crystal of 


NiO. 
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TUNNELLING FROM A MANY-PARTICLE POINT OF VIEW* 


J. Bardeen 
University of Illinois, Urbana, Illinois 
(Received December 16, 1960) 


Giaever' and more recently Nicol, Shapiro, and 
Smith? have observed the tunnelling current flow- 
ing between two metals separated by a thin oxide 
layer. The most interesting results are obtained 
when one or both of the metals are superconduc- 
ting, in which case they find direct evidence for 
a gap in the quasi-particle spectrum of the super- 
conductor. They were able to account for the 
data quantitatively on the assumption that the only 
relevant factor is the density of states in energy. 


This is to be expected if the transition probability 
for transfer of an electron from one side to the 
other is given by the familiar expression 
(27/R)|M "pp, where M is the matrix element 
and py the energy density of final states, and 

if it is further assumed that M can be treated as 
a constant. It is implied that M is not only in- 
dependent of energy for the small energy differ- 
ences involved, but is also unchanged when the 
metal goes from normal to superconducting. 
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However, it is not immediately obvious that these 
assumptions are justified. We give here a dis- 
cussion of tunnelling from a many-particle point 
of view and show that it is plausible to treat M 
as a constant in the interpretation of the experi- 
ments. 

We suppose that the barrier extends from x, 
to x», with metal a to the left of x, and metal b 
to the right of x,. Consider two many-particle 
states of the entire system, ¥, and ¥,,,,, which 
differ in the transfer of an electron from a to b. 
We suppose that ¥, and ¥,,,,, may be defined in 
terms of quasi-particle occupation numbers of 
metals a and b, so that W,,, differs from ¥, in 
the transfer of an electron from state m in a to 
state m in b, all other occupation numbers re- 
maining the same. Of course m may correspond 
to a normally occupied state, in which case there 
will be a hole in m in ¥,,». 

The quasi-particles do not correspond to plane 
waves, but to waves which are reflected at the 
barrier and which drop exponentially with dis- 
tance into the barrier region. For example, in 
the free-electron approximation for the normal 
state the wave function would be of the form 
(WKB approximation): 


-1/2 i 
>, 9. / ell Pyy *Pe)sin(p x +7), xSK (1a) 


- -1/2 i( pyy +p,2) 
Yn = 2C'0, | eo"? 


x 
xexo(-f 1b, I), % S*SK, (1b) 


Xa 


where C =(2p,. /L) is a normalization constant, 
and in the barrier region, |p, | =(2uU = Py” -p,"), 
where U(x) is the potential energy. We have taken 
units such that#=1. Beyond x,, we assume that 
the wave function representing the electron in 
state m drops smoothly to zero, instead of oscil- 
lating, so that it is not a good solution for x >x». 





mn 


- “lS an) J an*®o! , 


Thus we assume that W, is a solution of the 
Schrédinger equation with energy W, for x<x,, 
but there is a region to the right of x, where it is 
not a good solution. Similarly, we assume that 

V nn With energy W,,, is a solution for x >xq, 

but not for the region to the left of x, where the 
wave function for quasi-particle n drops to zero. 
Both ¥, and &,,, are good solutions in the barrier 
region x, <x <Xp. 

We form a time-dependent solution as a linear 
combination of ¥, and various final states, ¥,», 
by the usual method: 

-iWynt 


- -iWot 
Vv =a(t)¥ 0° +2 in mn mn? : 


(2) 


and substitute into the Schrédinger equation. 
This gives for the matrix element for the transi- 
tion 


* 
M nt f Wo" H-W Oe dr. (3) 


Since the integrand vanishes except over a region 
to the left of xg, we need to integrate only over 
this region. 

We may express the integral in a more sym- 
metric form by subtracting ¥,,,,,(H - W,)¥o> 
which vanishes to the left of x,. Since we are 
interested only in final states such that W,,,~ W,, 
the result may be written 


M J HY -¥ HW "ldr, (4) 


where the subscript a indicates that the integra- 


tion is to be taken over the region to the left of xj. 


The matrix element may be expressed in terms 
of that of the current density operator, J, in the 
barrier region as follows. We introduce into the 
integrand a step function S(x) which is equal to 
unity between a point x, to the left of the im- 
portant region of integration and a point x, in the 
barrier and which vanishes elsewhere. Then 
integrations with respect to y; and z; vanish and 
integration with respect to x; gives 


on 1 *— 2 2 * 
7 J fx, Ve ¥ mn * mn i %o S%)47 "°° Ay 


(5) 


where J,,,,,(x) is the matrix element of the x component of the current density operator defined by 


wy f...fte.* ’ * . ™ 
a ne oD) f f [¥,"ow 8%, -¥ OH 9" /ax,]Ox -x,)dr ++ -dr,. 
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According to our assumptions, J,,,,(x.)=0 and 
further, J),»(x,) is independent of position as 
long as x, is in the barrier region. Thus 
Mt os : “San y)- (7) 

It is easily verified that this method leads to the 
usual results for barrier penetration problems. 

The quasi-particle energy in a superconductor 
is E =(€?+A”)”, where e(k) is the normal state 
energy measured from the Fermi surface. In 
calculating the density of final states, it is simp- 
lest to take the conventions that f =1 and E is 
negative for the normally occupied states in the 
Fermi sea, so that for k<kp, the probability of 
a hole excitation is 1-f and the energy of a hole 
is -E. This is the procedure which has been used 
in the interpretation of the experimental data.'»? 
In each region of k space the density of states in 
energy in the superconductor differs from that in 
the normal metal by the factor 


p.=p, LIE! /(E? - a7)", (8) 


with p,=0 in the gap, |E|<A. Agreement with 
experiment is obtained if it is assumed that 
Mn (or Jmn) is the same for the corresponding 
transitions in normal and superconducting states. 
Usually coherence factors, which have a marked 
effect on transition probabilities, are introduced 
in calculations of matrix elements between quasi- 
particle states of a superconductor.* To see why 
such factors are not expected to occur in tunnel- 
ling, one must consider the pairing in the vicinity 
of the barrier. Coherence factors would be in- 
troduced if one simply paired complex conjugate 
wave functions of the type (la), (1b), so that 
pairing extends into the barrier region. How- 
ever, if one looks at the problem more closely, 
from the viewpoint of the more general Gor’kov 
equations‘ which allow for a variation of the energy 
gap parameter with position, one sees that A will 
drop to zero very rapidly in the barrier. In effect 
electrons in this region are not paired and the 
wave function is essentially the same as in the 
normal state. 





For an effective interaction v(f,,f,), the posi- 
tion-dependent energy gap function is defined by 


A(f,, f,) - (N -2| v(F, )w(F,) IN) av v(F,, B,). (9) 


Simply from the fact that y is very small in the 
barrier region, one expects that A(T,,f,) is small 
when either f, or F, is in the barrier. Note that 
the tail in the barrier of the wave function of a 
typical electron at the Fermi surface is very 
much smaller than that of one of the few elec- 
trons moving normal to the interface which has 
an appreciable probability of penetrating. How- 
ever, A(f,,F,) is expected to rise very rapidly to 
normal values in the superconductor, A(?,,f,) 

~ 6(F,-F,)A, with A =const. 

Since, according to (6), the matrix element de- 
pends only on the wave function in the barrier, 
it would be expected to be the same as for the 
corresponding transition in the normal state. 
Actually, the wave function in the barrier would 
be changed slightly because of the difference in 
quasi-particle energies in normal and supercon- 
ducting states, but this would have a negligible 
effect on the matrix elements. Thus the only 
significant difference in the tunnelling current 
comes from the density of states factor. 

The method described here can also be used 
for calculating the tunnelling current between the 
valence and conduction bands of a semiconductor, 
as observed in the Esaki diode. It can be gen- 
eralized to take indirect transitions into account. 

The author is indebted to Dr. W. A. Harrison 
and Dr. J. R. Schrieffer for discussions of the 
theory of tunnelling in a superconductor. 
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ELECTRIC QUADRUPOLE SPLITTING AND THE NUCLEAR 
VOLUME EFFECT IN THE IONS OF Fe*’t 


S. DeBenedetti, G. Lang, and R. Ingalls* 


Carnegie Institute of Technology, Pittsburgh, Pennsylvania 


(Received December 19, 1960) 


A study of the Mossbauer spectra of Fe*” in 
chemical compounds has revealed some striking 
regularities. Typical transmission curves for 
Fet* and Fet** ionic compounds relative to an 
unsplit room-temperature stainless steel source 
are shown in Figs. 1 and 2, respectively. Both 
patterns exhibit a shift of the center of gravity 
from the zero of velocity, and the curve of Fet* 
shows the two absorption lines characteristic of 
a pure nuclear electric quadrupole interaction. 
The quantities 6 and AE, which characterize 
these effects, are defined in the figures; their 
values for various compounds are shown in 
Table I. It is noteworthy that the value of 5 and 
the low-temperature value of AE are, for a given 
iron ion, relatively independent of the chemical 
combination. 

The center of gravity displacement, or chemi- 
cal shift,! arises from two mechanisms: (a) the 
second order Doppler shift,?** which is caused by 
lattice vibrations, and is a function A(7T5,7q,45,9q) 
of the source and absorber temperatures and 
Debye temperatures, and (b) the nuclear volume 
effect,*»®> which represents the change, from 
source to absorber, in the modification of the 
nuclear transition energy caused by electrons 
overlapping the finite nucleus. The total shift in 
a M&ssbauer spectrum may be written: 


2p 2p 2 2 
6 «A(T, TY Oy a) +B(R, “8, 1¥, (0) - _ (0)], 
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FIG. 1. Mossbauer spectrum of the Fe** ion in an 
absorber of FeSo,°7H,O at nitrogen temperature using 
a room-temperature stainless steel source. The pat- 
tern exhibits the chemical shift 6 and the electric quad- 
rupole splitting AE of the excited state of the Fe*’ nu- 
cleus. The velocity is positive for the source approach- 
ing the absorber. 
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where Re and Rg are the charge radii of the nu- 
clear excited and ground states. ¥2(0) is the 
s-electron density at the position of the nucleus. 

B is a constant which is 1.76 x10° cm/sec for 

iron when R2P and ¥2(0) are expressed in atomic 
units. p=(1-a?Z*)” is a relativistic factor® equal 
to 0.982 for iron. 

An estimate of Rg/Re may be made if one neg- 
lects the influence of the crystalline surroundings 
and uses wave functions® for the free iron ions. 
Using (1) and the room-temperature shifts (which 
are less sensitive to differences in 6 than are the 
low-temperature shifts), we obtain 


2p 2p 2 2 
=B(R , “a, Ie, atte (0)}. (2) 


—, - 


6 -6 
++ 0 +44 


We have neglected the small difference A,,-A,,,. 
Watson’s wave functions yield %442(0)< ¥,44°(0); 
the extra electron in Fe** apparently shields the 
3s electrons slightly. Equation (2) yields Rg/Re 
=1.001. 

Watson®,” has also calculated wave functions 
for the outer electron configurations 3d°4s? and 
3d* for iron. Using an equation similar to (2) one 
obtains the shifts expected for each case: 


Ore =-0.178 cm/sec for 3d*4s?, 


=+0.182 q 
ore +0.182 cm/sec for 3d® 
When metallic iron is used as an absorber with 
the stainless steel source at room temperature, 

a shift of +0.01 cm/sec results. Since this value 

is midway between the above calculated values, 
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FIG. 2. Mossbauer spectrum of the Fe+++ ion in a 


FeCl; absorber at nitrogen temperature using a room- 
temperature stainless steel source. The velocity is 
positive for the source approaching the absorber. 
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Table I. 


A list of quadrupole splittings AE and chemical shifts 5 relative to the emission line from a stainless 
steel source at room temperature. 











my 








Absorber Temperature AE (cm/sec) 6(cm/sec) References 
Fe,0; Room 0.024 + 0.003 0.047 + 0.003 a 
Fe,0; Room . 024 0.050 b 
Fe,(SO,)s Room 0 0.064 b 
Fe,(SO,)3 Room 0 0.055 + 0. 005 this work 
Fe,(SO,4)3 Nitrogen 0 0.065 + 0.005 this work 
FeCl, Room 0 0.045 + 0. 005 this work 
FeCl, Nitrogen 0 0.065 + 0.005 this work 
FeCl,°4H,O Room 0.300 + 0.005 0.135 + 0.005 this work 
FeCl,°4H,O Nitrogen 0.310 + 0.005 0.145 + 0. 005 this work 
FeSO,°7H,O Room 0.320 + 0.005 0.140 + 0.005 this work 
FeSO,°7H,O Nitrogen 0.360 + 0.005 0.150 + 0.005 this work 
Fe(NH,) (SO,)2°6H,O Room 0.175 + 0.005 0.140 + 0. 005 this work 
Fe(NH,),(SO,).°6H,O Nitrogen 0.270 + 0.005 0.150 + 0.005 this work 
FeF, Room 0.268 c 





See reference 4. 


by, Solomon, Proceedings of the Allerton Park Conference on Mossbauer Effect, University of Illinois, 
Urbana, Illinois, June 5-7, 1960 (unpublished), Sec. IIc. 


°G. K. Wertheim, Phys. Rev. 121, 63 (1961). 


it strongly suggests that metallic iron has an 
outer electron configuration equivalent to 3d74s'. 
Although the electric quadrupole splitting seems 
to be temperature dependent, it is roughly the 
same at low temperatures for all the Fe** com- 
pounds listed in Table I. This suggests that the 
relevant field gradient is mainly caused by the 
ion alone. The free Fe** ion has an outer elec- 
tron configuration of 3d° and is in a °D state. 
Five of these 3d electrons have their spins aligned 
and together form a spherically symmetric dis- 
tribution of charge. The sixth electron has oppo- 
site spin and also a choice of the five possible 
values of mj. The crystal symmetry dictates the 
m , combination the ion will actually take. For 
sufficiently low symmetry the lowest energy state 
will be a linear combination of the two states with 
m, equal to 2 and -2. Assuming this, and neg- 
lecting all other possible effects of the crystalline 
environment, it is possible to calculate the ionic 
field gradient at the nucleus to within the factor 


(1-y), where y is the Sternheimer antishielding 
factor. Taking AE to be 0.3 cm/sec one obtains 
|Q(1-y)!=0.1 b, where Q is the nuclear electric 
quadrupole moment for the excited state of Fe*”. 

The temperature dependence of the quadrupole 
splitting could be caused by thermal excitation of 
the Fe** ion into a state which yields a lower field 
gradient at the nucleus. For sufficiently rapid 
ionic transitions, the nucleus will respond to the 
time average of the field gradient. This interpre- 
tation indicates that the level splitting of the Fe** 
ion in Fe(NH,),(SO,),-6H,O is sufficiently small 
that significant excitation occurs at room tempera- 
ture. 

It is noted that there is very little quadrupole 
splitting in the case of Fe+t*+. Presumably there 
is negligible field gradient from this ion, since 
its 3d° electron configuration gives a °S state. 

The splitting that does occur‘ is therefore caused 
by field gradients from the neighboring ions rather 
than from the Fe*** ion itself. 
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DETERMINATION OF THE PARITY CHANGE IN INELASTIC SCATTERING 
FROM THE ANGULAR DISTRIBUTION* 


A. J. Kromminga and I. E. McCarthy 
Department of Mathematical Physics, University of Adelaide, Adelaide, South Australia 
(Received November 21, 1960) 


Considerable success in determining the angu- 
lar momentum transfer in direct-interaction 
inelastic scattering has been achieved by com- 
paring the angular distribution in the approximate 
range 20° to 90° with |), (KR)? or the improved 
theory due to Austern, Butler, and McManus.’ 
Further information that one would like to obtain 
is the parity change of the nuclear state. 

The distorted-wave Born approximation has 
been shown to give rather detailed fits to experi- 
mental data whether the final state is a single- 
particle’ or collective® state. In a simplified 
distorted-wave Born approximation, Glendenning* 
found theoretical examples of angular distribu- 
tions in which the forward differential cross sec- 
tion was very small for odd-parity change and 
very large for even-parity change. He suggested 
that this could be a general result which would 
make it possible to determine the parity change 
by a very simple experiment. 

It is shown here in the distorted-wave Born 
approximation that, for reasonably small mo- 
mentum transfers, the forward differential cross 
section is always small for odd-parity change 
and large for even-parity change. This result 
provides a very powerful tool for nuclear spec- 
troscopy. 

The matrix element for inelastic scattering in 
the distorted-wave Born approximation may be 
written 


mer f dr oe M({k, r)o*(k’, rv (r), (1) 


where y“ and gy‘ are the wave functions of the 
incident and scattered particles whose momenta 
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are hk and hk’, respectively; ¥ (Tf) contains all 
the properties of the initial and final nuclear 
states and the interaction. In the case of single- 
particle excitations, 


V@)=f ary, EWE FF), (2) 


where y; and y are the initial and final single- 
particle wave functions. 

In the case of rotational excitations by spin- 
less particles, in the notation of reference 3, 


Wr) wp l¥, "(0,,69)1v) ¥, "ay au/ar), (8) 


where 6,,6, are the orientation angles of the nu- 
cleus with respect to the space-fixed frame; 2 
is the angular part of r; vg and v; are the final 
and initial wave functions, respectively; and L 
is the orbital angular momentum change in the 
reaction. 

The point is that the parity of ¥(F) in all cases, 
whether the incident and scattered particles have 
spin or not, is determined by the parities of the 
initial and final nuclear states; that is the parity 
change of the nucleus. 

The remaining factor in ¥(F) has even parity 
including the pairs of spinor wave functions and 
spin- dependent interactions in the case where 
the scattered particles have spin. 

The functions g‘ and g‘~* may be written 


g(r) =X (al+ vil exp(io,) R (er P, (cosé), 


gy *(r) = 2, (21 +1)(- iy! exp(io )R'r)P (cosé’), (4) 
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where @ and 6’ are measured from the direction 
of k and k’, respectively. 

Assuming that k and k’ are not very different 
and that the optical model potentials seen by the 
incident and scattered particles are similar, we 
notice that, if the directions of k and k’ are the 
same (forward scattering), 


ge) = g(r, 7- 4), 
gy *(r) = g(r, 7- 6). 


Hence the product g(r) y*(F) has approxi- 
mately even parity. 
If ¥ has odd parity (nuclear state undergoes 
parity change), 
om (k’ =k) = 0. 


If ¥ has even parity (nuclear state undergoes no 
parity change), 
Mt (k’ 2 k) #0. 


When ¢‘ and y‘™ are approximated by plane 
waves, the matrix element for forward scatter- 
ing becomes 


- _ ” 1 
om (k’ > k)x f vdr f exp(iKr cos@)d(cosé) 


0 0 


2n 
x f ¥@ag, (5) 


where K=k -k’. If the contribution from L =0 to 
J¥ (F)dg is small, then m(k =k) =0. 

When optical model wave functions are used 
for y and g**, there is a strong 6 dependence 
of the modulus of yg“ and y“*. The integrand 
cannot be written in the form (5), since there is 
another 6-dependent factor. Hence if ¥ has even 
parity, 9on(k~=k) will not in general be small. 

The @ dependence may be understood by con- 
sidering yg. Let us write 


g(r) =A, (r, 6) exp[iS(r, 8)]. (6) 


Then Ay decreases slowly as 6 decreases from 

7, and, for nucleons of energy less than about 40 
Mev, increases sharply as 6 goes from about 

15° to 0° (the focus). It has been shown’ that the 
slow decrease from @=17 may be approximately 
absorbed into the factor exp(iKr cos) by making 
Ky complex. However, the focal contribution can- 
not be absorbed in this way and is responsible 

for the large value of om (k’ =k) when W has even 





parity. 

Glendenning’s example suggests that, for nu- 
cleons of about 20 Mev, it would be possible to 
distinguish rising from falling differential cross 
sections by measurements down to 10°. 

A crude estimate of the condition for the for- 
ward differential cross section to be small can 
be obtained from the approximation of reference 
5 which has proved very successful in describing 
angular distributions in inelastic scattering of 
@ particles: 


y(®) =A, (r, 8) exp(ikr cos), 
gy *(r) =A, ,(r, 1-0) exp(-ik’ cosé), 


ot (k’ = k)« ('d(cosé) A, (79 BA, , (7597-8) 


x exp(iKr) cosé) f dy v(F,). (7) 


Here we have restricted the radial integration 
to r=7o. 

The non-nuclear part of (7) is expected to have 
even parity if K is sufficiently small for the sine 
part of exp(iKr, cos@) to be always small com- 
pared to the cosine part. sinKr,<0.2 for r,=5f 
if K<0.04 f-*. For 20-Mev neutrons, this cor- 
responds to |Q|<1.6 Mev. For higher energies 
1Q| may be larger (1Q| nay «E™*), This condi- 
tion is easily satisfied by low-lying states of 
medium nuclei. Exact calculations should give 
a better idea of the range of applicability of the 
rule. 

The result should apply also to reactions where 
the final and initial nuclei have the same mass 
such as (n,p) reactions, since the Coulomb po- 
tential does not severely modify the optical 
model wave functions in the interaction region. 
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ROLE OF THE PERIPHERAL INTERACTION IN PROTON-PROTON COLLISIONS 


Franco Selleri 
CERN, Geneva, Switzerland 
(Received December 19, 1960) 


It has been suggested in a recent paper by Bon- 
signori and the writer’ that the inelastic scattering 
of any particle on a target nucleon may be dom- 
inated by the exchange of one single virtual pion. 
Analyzing some known experiments,”»* it was there 
shown that in the Bev region there is good evi- 
dence for such a mechanism, when the incoming 
particle is a 7 meson or a nucleon. Further ex- 
periments* have shown the same qualitative evi- 
dence which consists essentially in a low-energy 
peaking of the laboratory spectrum of the recoil 
nucleon. It seemed therefore desirable to go 
further and calculate the cross section for the 
process 


(1) 


p+p—p+nen* 


under the assumption of dominating peripheral 
interaction. We wish to stress that this dominance 
is assumed throughout the whole physical region 
and not only for very small momentum transfers 
as has been done in recent work.® The reason 

for choosing reaction (1) is that final-state inter- 
actions between the nucleons are negligible at 

high energy, while the (3,3) interaction of the 

pion is likely to involve only one of the nucleons 
and thus not to destroy the typical effects of the 
one-pion exchange. A typical peripheral diagram 
for this process is shown in Fig. 1. p, and p, 

are the 4-momenta of the incoming protons; q, q,, 
and g, those of the outgoing pion, proton, and neu- 
tron, respectively. Together with this, one should 
consider the other three diagrams obtained by ex- 
changing the initial protons and the final nucleons. 





m 
20 E EEE 
( q~P, b2 G a.) 





Si = -(27)” 


where 


U(qp)y54( Pe) ” 
M( P,P) = (a, bJro wal -A + iqgB \u( p,)- (3) 


In (2), wg(E 7) is the zero component of the four- 
vector q(l), Pf =9+4,+42, P;=P,+p2, and V2G, 
is the rationalized and renormalized coupling 
constant for emission of a charged pion. The 
factor 2~“? and the minus sign are typical effects 
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FIG. 1. Peripheral diagram contributing to proc- 
ess (1). Full and dashed lines represent nucleons 
and pions, respectively. 
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However, simple isotopic spin arguments show 
that the separate contribution to the total cross 
section for process (1) of the two diagrams in 
which the 7° and the neutron emerge from the 
same vertex is one ninth that of the diagram of 
Fig. 1. Therefore we can consider them of sec- 
ondary importance and neglect them in our cal- 
culation. Furthermore, we shall be forced to 
calculate the nucleon-nucleon-pion vertex and 
the pion propagator to lowest order, and to sub- 
stitute for the four-particle matrix element its 
expression known phenomenologically for the 
case in which all four particles are on the mass 
shell. We will see a posteriori upon comparison 
with experiments that these approximations prob- 
ably are not too bad. 

In the stated approximation the S-matrix ele- 
ment for process (1) can be written 


(2) 





of the identity of the initial protons.® In (3), A and 
B are the well-known invariant amplitudes for 
pion-nucleon scattering (in our case in the state 

of isospin 3/2) as defined for example by Chew 

et al.,” and m (yu) is the nucleon (pion) mass. From 
(1) and (3) it is in principle possible to calculate 
any cross section at any energy. We choose to 
calculate do/dT, T being the lab kinetic energy of 
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the neutron. The importance of this variable was the Birmingham group.” At this energy, the 
stressed in reference 1 and is essentially due relative energies of the 1 and the proton are in 
to its simple relation with the 4-momentum trans-___ the region where the (3,3) resonance dominates. 
fer to the target nucleon: (p,-q,)?=2mT. In We have therefore kept only the J=3/2 P-wave 
’ performing the actual calculation we had in mind part of A and B. For future reference we write 
a comparison with the 970-Mev experiment of down the result of our calculation: 
a. 8 f £7), oR (w)o(w)[a(T) +0(T,, w) +c(T, w) dw, (4) 
os i M mM+U 








where b* = 1.7(167m”)"'G,? = 0.08, p; is the lab 


2 2 
momentum of the incoming proton, w =[-(q+9,)"]” In the preceding formulas A,° and A,° are defined 


is the total energy of 7* and p in their c.m. sys- by 

tem, and o(w) is the 7*p cross section at a c.m. A,’ =(p, -9,)? =m? + 2m(T; -T) -w’, 
energy w.® R, a,b, and c are the following func- 

tions: A,’ =(b,- 42) =2mT, 


R(w) = 4[w* - 2w?(m? + p?) +(m? - u?)?] 2, 5 T; being the lab kinetic energy of the incoming 
(w) = al ( we) + wy] (5) protonandE,’, E,', 6; the energies of the 


a Snr 


A,’ initial rotons and the le between their direc- 
@(T) =-—s*—ay (6) P ang 
(7) (A,* +p*)°’ tions, all in the c.m. system of the 7* and the 


A.? final proton. These three quantities can be ex- 
b(T, w) a eT , (7) pressed as a function of W, total c.m. energy for 
A > 
ae reaction (1), w?, 4,?, and 4,” through the follow- 
(8) ing relations: 





: ¢ 
OCP) = TR MA)” 


E, "= (w? +m? +4,°)/2w , (10) 
1 
o= (w mS,” +e, * +m)]|? cos6 ; E ‘= (w? +m? +4,°)/2w ’ (11) 
2 
: 29. - = 
| +(w+m)[(E, ’-m)(E, '-m)]|” siataees, . c0s6; =(P,*P2)/P:b2, (12) 
Pe A . (9) where p, and p, are three- momenta correspond- 


PERIPHERAL MODEL 


STATISTICAL 
THEORY 


L_— 





tp © 





FIG. 2. Experimental 
energy spectrum in the 
lab system for neutrons 
from reaction (1) as com- 
pared with statistical 
theory (dashed line) and 
our formula (4) (full line). 
































VoLUME 6, NUMBER 2 


PHYSICAL REVIEW LETTERS 








JaNuary 15, 1961 





ing to Ep,’ and Ep,’ respectively, and 


P (2a - A,”)(2a@ - A,”) 


P,P, =-@ 4w2 ’ (13) 





with a=W?/2-m?. Finally we have for the upper 
limit of integration 


7? 27\2]" 
et) [rt we Zowene -amaye( +) 
(14) 


The predictions of formula (4) at 970 Mev are 
shown in Fig. 2 and compared with the experi- 
mental spectrum.® As is seen, the agreement 
is good while a statistical spectrum fails to ex- 
plain the three main features of the experimental 
distribution; these are: the strong peaking at 
low (~50 Mev) energy, the broad bump at high 
(~500 Mev) energy, and the minimum between 
them. These features are, however, quantita- 
tively predicted from our formula (4) and are 
typical effects of the peripheral interaction (see 
reference 1 for a more detailed discussion). 

Theoretically, the total cross section’® is pre- 
dicted to be 12 mb; its experimental value is 
16.4+0.7 mb.” This indicates that something 
more than the diagrams of Fig. 1 contributes to 
process (1). However, this “something,” what- 
ever it is, is not more than 30% of the peripheral 
contribution. In view of these results it seems 
very desirable to have a deeper comparison be- 
tween this model and the experimental data." 
The author is indebted to Professor A. Stang- 












hellini and Dr. E. Ferrari for very useful sug- 
gestions, and also to many theoreticians of CERN 


for discussions. 
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INTERMEDIATE BOSON PRODUCTION IN PION-PROTON COLLISIONS* 


Norman Dombey 


California Institute of Technology, Pasadena, California 
(Received December 7, 1960) 


The current-current hypothesis has been ap- 
plied in recent years to the study of the weak 
interactions. The conserved vector current’ hy- 


pothesis and the pionic character of the divergence 


of the axial vector current’ have evolved in its 
wake. So has the idea of a boson field which me- 
diates all weak interactions.* This note deals only 
with the vector part of the current, leaving for 
another time similar calculations with the axial 
vector part; and examines the possibility of pro- 
ducing the boson W, in pion-nucleon collisions, in 
case of W having relatively low mass.* 
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For definiteness let us consider the reaction 


1 +p~-W* +p, (I) 


and concentrate solely on the vector part of the 
interaction. Then, on inverting incoming and 
outgoing particles, it is apparent that the process 
is almost exactly the same as electroproduction 
of pions (i.e., production of pions by photons off 
the mass shell). The differences are: (i) The 
coupling constant is not the electromagnetic coup- 
ling constant. (ii) The four-momentum of the 
“photon” is timelike, not spacelike. But the anal- 
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ysis of electroproduction in papers A and B* is in 
fact relevant provided that instead of e? = 1/137, 
we write’ e? =(1/nV2 )(N?/M?)x107°, where N is 
the boson mass and M is the proton mass. This 
simple substitution is the content of the conserved 
vector current hypothesis. The problems raised 
by (ii) imply the extrapolation of the analysis in 
A and B to a new region. 

The total amplitude 9% may be written in the 
form 





mm =AM , +BM_+CM_.+++++FM_, (1) 


A B Cc F 
where the six functions A, B, «++, F depend on 
the invariants 
_ Pek 7 q:k 

v="? "RB" OM’ 
here P=3(p,+ 2); p;, p2 are the final and initial 
momenta of the proton; gq, k the momenta of the 7 
and W. M4, °::, Mp are defined by Eq. (4) inA 
in terms of y matrices. It is also possible to 
write sm in terms of two-component Pauli spin 
matrices: 


k? = -N?; (2) 


6 
m=F=2,2.5., (3) 

ae ee 

i=1 
where the ©; and the linear connection between 
the $; andthe A, B, ---, F are given by Eqs. (8) 
and (9) in B. The cross section in the center-of- 
mass system is then given by 


do \kl M? 


an Tar we Al slDl’, (4) 


where W is the total energy in the c.m. system, 
and the matrix element includes the sum over 
final and average over initial spin states. 

As the energy for reaction (I) is much above the 
3-3 resonance of the pion-nucleon system, the 
calculation first will be done using only the Born 
terms. Here 


A=-ef Pw ni(, Sretvn :) 

















BY "Be 
- if, vrai), 
Can, ace 7 ie) ' 
| D=u, FPR) (- v* * ) ; 





E=F =0, f?=0.08, #,’= 1.786 /AM), (5) 





F,? and F,? are the Hofstadter form factors of 
the proton. Note that the form factors are of 
opposite sign in argument to those determined 
experimentally. The result of the calculations 
done with F,P(k?) =F ,P(k*) =1 for varying N at 
fixed energy is shown in Fig. 1. The cross sec- 
tion varies slowly with the energy for a given W 
mass. 

It has been pointed out by Frazer and Fulco® 
that a 7-7 resonance in the state J=1, J=1 would 
give theoretical agreement with the experiments 
for the isotopic vector form factors. A more 
recent paper by Bowcock, Cottingham, and Lurié’ 
fits the position and width of such a resonance to 
the data on 7-N scattering and electromagnetic 
nucleon structure. 

From the dispersion relations (with subtractions 
where necessary) 


ipo 8) (eat! 


V 
F, Os). Pa 


_ , (é=1 or 2) (6) 
1 4p? 
with 


¥. \- , V , 
g, (t= IF (Mle, (tN 5 (7) 


where the [g;(t’)], are given in the papers on 
nucleon structure neglecting the 7-7 resonance,® 
the real and imaginary parts of F;"(t=+N*) can 
be calculated. The pionic form factor is given 
(for ¢ not too large) by 


t 2 
Y 
IF (e)I" = (f-t P+y7v*? v=it-n', (8) 


when the m-7 scattering amplitude is written in 
resonance form (for J=1, J=1) 


faa 
11 +, oe 


lor 
8e 

x 6 } 
(lo-*Fem*) 4 | 
ob 





i 1 — 1 l 


fe) 
354 5 6 7 8 9 10 
N (in units of 2) 


FIG. 1. Cross section o as a function of boson mass 





N in uncorrected Born approximation. Momentum of 
pion in c.m. system is 980 Mev/c. 
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FIG. 2. Pionic form factor about resonance region. 


Curve I given by Frazer and Fulco. Curve II given by 


Bowcock, Cottingham, and Lurié. 


In reference 7, the values for y and tare given 
as 


y =0.376y"?, t =22.4p?. (9) 
The two alternative curves for |F,(¢) |? are given 
in Fig. 2. The multiplicative factor (taking the 
charge and magnetic moment form factors to vary 
similarly in the low-energy unphysical region as 
well as the physical) to the curve in Fig. 1 is 
given in Table I. 

An estimate of the effect of the 3-3 resonance 
was attempted using the dispersion relations for 
A, +++, F given in A and B and replacing the reso- 
nance by a delta function, i.e., 


(1/7)ImA (x) =A 4q(x)6(1 - w/w), etc., (10) 
where x=v-vg, w=W-M, w, refers to the posi- 
tion of the resonance, and A,,(x) is the 3-3 pro- 
jection of the Born term for A, and is given in 
paper B. As the 3-3 resonance occurs well into 
the unphysical region for reaction (I), the ampli- 
tudes A,,, ---, F,, must be analytically continued. 
If this is done and the calculation carried out for 
the total cross section, a very large increase is 
obtained (>100). There is no reason, however, 

to take this calculation seriously; but it should 
be viewed as an indication that the cross section 
may yet be larger than that given by the Born 
term with or without corrections. 
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Table I. Multiplicative factor given by 7-1 reso- 
nance to curve in Fig. 1 for various values of N.* 








N 3.5 4 4.5 5 6 7 9 
Frazer and 

Fulco 5 —{ + _ J <1 <1 
Bowcock, 

Cottingham, 

and Lurié 20 45 150 60 15 § ~1 





4° (2) is neglected. 


So, in summing up, it can be said that if the 
Bowcock, Cottingham, and Lurié fit to the 7-7 
resonance is correct, then the cross section 
should be of the order of 10~*° cm? for values of 
the W mass between four and five pion masses. 
It is also possible, if the 3-3 resonance is taken 
into account, that the cross section is greater 
still; and it must be remembered that the axial 
vector contribution has been completely neglected. 

I would like to thank Dr. Leon Landovitz for 
suggesting this problem, Dr. G. C. Wick for his 
hospitality at Brookhaven National Laboratory, 
and Dr. Murray Gell-Mann and Dr. Fredrik 
Zachariasen for comments and advice. 
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BINDING ENERGY OF THE TRITON 


G. Derrick 
University of Manitoba, Winnipeg, Canada 


and 


D. Mustard and J. M. Blatt, 
New York University, New York, New York and the University of New South Wales, Sydney, Australia 
(Received December 9, 1960) 


We have made a variational calculation of the 
binding energy of the triton, using the same two- 
body nuclear force law which was used by Brueck- 
ner and Gammel in their calculation of the pro- 
perties of nuclear matter,’ and assuming, as 
they did, that there are no specifically many- 
body forces. Our results disagree with the ex- 
perimental data on the binding energy of the triton 
and on the Coulomb radius of He*, thereby ex- 
cluding the Brueckner-Gammel force law. 

Of the 10 different states which contribute to 
the ground-state wave function of the triton,” we 
have included 8. The two states omitted have 
fully antisymmetric “internal” wave functions 
(states Nos. 2 and 5 of reference 2), and thus 
high kinetic energy. 

The trial wave functions were chosen as fol- 
lows. For each state i (i=1, 2, ..., 10), let 


y,= expl-y,(r-75)], (1) 


where 7, is the core radius, and y; is an adjust- 
able parameter. Let aj3, aj5, a;7 be further 
parameters, and let 


u(r)== Ly +a,9(y )+a.0(9,)° 
+a-(y)" -(1 +a. +a,,+4-)(y)°} (2) 


Then the trial wave function for the dominant, 
symmetric, *S,, state has the form 


Yy =A yt (7 yp) (Fog), (751) +A Uy (7 yQ)ley (75), (75), 
(3) 
where A,, A, are parameters, and u,(r) differs 
from u,(r) only through the parameter values 
Yir Figs 415, F17- 
Wave functions of “mixed” symmetry, e.g., 


the mixed-symmetry *S,, state (state 3 of refer- 
ence 2) are constructed as follows: Let 


z= exp[-26 (r - 10) 1, (4) 


where 6; is a parameter, let aj4, @j6, and a;g 
be additional parameters, and let 





1 
vir) = 12 +444 (2 Pra, els 


+a. 


_- 
+a,,(z )*- (l+a,, +a, 


5 
Then the internal part of the final wave function 
for the mixed-symmetry *S,,. state has the two 
components 


3,1 =Ag(3)~ M2 Luis (7 2) [U5(723)M5 (751) + 45 (725)Y5 (51) ] 
~ 20, (7,2)ug(%25)45(75,)}, 
3,2 = Ags (712) [v5 (77p5)45(7g1)-4g(725)¥ (751) ]- (9) 


Trial functions for P states and D states are con- 
structed in an analogous fashion, with the appro- 
priate factors in front so as to give the correct 
dependence on Euler angles and spins and the 
correct asymptotic behavior in the singular con- 
figurations (equilateral triangle, and triangle 
degenerating into a straight line).* However, the 
doubling-up of wave functions for one state, as 

in (3), is employed only for the dominant *S,, 
state, not for any other states. 

The matrix elements of the Hamiltonian, using 
these states, involve sums over spin and isospin 
indices, integrals over Euler angles, and inte- 
grals over 7,,, %23, %3,- The spin sums and Euler 
angle integrations are done analytically,** the 
integrations over 7,,, 73, %3, numerically.* De- 
pending upon the fineness of the net and upon the 
number of states used, the numerical integrations 
take between 10 minutes and 1 hour on the IBM 
704. The number of netpoints ranges between 
9000 and 30000, and the upper limit on 7; between 
R=6.4x10~ cm and R=9.6x10-'% cm. The domi- 
nant *S,, state integrals are done with a much 
finer net which gives 5-figure accuracy; the other 
integrals have numerical errors of the order of 
1%. The numerical error in the expectation value 
of the Hamiltonian is estimated to be less than 
0.3 Mev. 

After varying some 40 parameters, we find for 
the ground-state energy of the triton 


E,<- 10 Mev, (7) 
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and for the Coulomb radius 
R. =3.2x1075 cm. (8) 


These values disagree with the experimental 
values, -8.49 Mev and 2.26x10-'* cm, respec- 
tively. Since this is a conventional variation cal- 
culation, any improvement in the wave function 
would increase the disagreement in the energy. 

We do not believe that the upper bound (7) for 
E, is close to the mathematically correct ground- 
state energy. Every time we varied new para- 
meters, or revaried ones varied much earlier, 
the bound on E, decreased another few tenths of 
an Mev. We merely quit varying parameters 
when the discrepancy with experiment became 
sufficiently large to be significant. 

The D-state probability is about 7%; no direct 
comparison with experiment is possible for this 
quantity. 

The P states are completely unimportant. 
States 4, 6, and 7 of reference 2 have probabi- 
lities 3.8x107*, 0.8x10-*, and 1.9x107*, re- 
spectively, adding up to a total P-state proba- 
bility of 6.5x10-*. The contribution of the P 
states to the binding energy is less than 100 kev. 

The L-5 force contributes significantly, but 
not through the P states. Rather, the matrix 
elements of the L-S force between the D states 
must be included in the calculation, and these 
matrix elements contribute several Mev, ina 
repulsive direction. 

The sharp variation of the trial functions near 
the core radius, especially for state 8 of refer- 
ence 2, requires a reasonably fine net for the 
numerical quadratures. 

The large number of states required, and the 
sensitivity of the result to the values of many of 


the parameters in each state, makes it appear 
doubtful that any simple approximation can give 
accurate results in a three-body problem with 
tensor forces. By analogy, simple approxima - 
tions in the many-body problem may well turn 
out to be less accurate than supposed. 

We conclude that the combination of the Gammel- 
Brueckner potential with the assumption of super- 
position of 2-body forces, is inconsistent with 
experimental data on the triton and He*. This 
does not necessarily contradict the result of re- 
ference 1, since one and the same potential may 
disagree with H* and He® and yet agree with nu- 
clear matter. However, it is not impossible that 
the accuracy of the Brueckner theory of nuclear 
matter is somewhat less good than had been hoped. 

We are very grateful to the U. S. Atomic En- 
ergy Commission and the Institute of Mathemati- 
cal Sciences, New York University, for travel 
allowances and a very generous allocation of ma- 
chine time, which alone made this calculation 
possible. We greatly appreciate many helpful 
and informative discussions with Professor R. D. 
Richtmyer, Professor E. Isaacson, Professor 
Max Goldstein, and Mrs. B. Weitzenhoffer, and 
the unstinting and generous help given to us by 
them and by all the members of the staff of the 
Atomic Energy Commission Computing Center 
at New York University. 
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THE REACTION y +y —v +7! 


Murray Gell-Mann 
California Institute of Technology, Pasadena, California 
(Received December 21, 1960) 


In a current-current picture of the weak inter- 
actions, a coupling term of the form (ev)(ev) 
arises.'»? The contribution from charge-exchange 
currents is an interaction Lagrangian 


ler, +y,)v ley (1 +7_)v] + ; (1) 
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in the case of local coupling. [If there are charge- 
retention currents involving leptons, other terms 
may occur, of course, conceivably even cancel- 
ling (1).] 

In reference 2, the cross section for neutrino- 
electron scattering is calculated using the inter- 
action (1); the result is too small for an experi- 
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mental test to have been possible up to the present 
time. 

Pontecorvo® has suggested that if the coupling 
(1) exists, it may be of astrophysical importance 
because it would induce reactions leading to en- 
ergy loss in stars. Recently* consideration has 
been given to processes like 


e-+e* + +7, (2) 


which is similar to electron-neutrino scattering, 
and 


Y+y eve. (3) 


It is interesting that the latter process is not, 
in fact, induced by the local coupling (1). The 
reasoning is as follows: 

By a Fierz transformation, we may express the 
coupling as 


Sfrgl+relrg(t+yv], (4) 


and we must then calculate the matrix element of 
2y qe and of 2ygys5e between the vacuum and a 
state of two photons. The first vanishes by charge 
conjugation invariance. The second is analogous 
to the amplitude for the decay of a pseudovector 
or a pseudoscalar meson into two photons. Since 
a spin 1 meson cannot decay° into two photons, 
only that part of the amplitude corresponding to 
a pseudoscalar meson survives. In the rest 
frame, only the time component can contribute. 
Thus the entire matrix element of 2y)(1+y5)e is 
equal to P, (the total four-momentum) times a 
pseudoscalar bilinear in the field strengths of the 
two photons. But the Dirac equation for the 
neutrinos tells us that DPyyg(1+y5)v gives zero. 
Thus local weak interactions, to lowest order 
in G, give a vanishing rate® for reaction (3). 


Our proof does not hold, however, for the case 

of an intermediate charged vector boson? medi- 

ating the charge-exchange current-current inter- 

actions. That case must be investigated further. 
The process* y+y -v+i+y is not forbidden 

for local coupling, although the contribution from 

@yqyv5¢ Vanishes by charge conjugation symmetry. 

The remaining amplitude, coming from 2yge, 

can be calculated from the amplitude’ for the 

corresponding process where the neutrino pair 

is replaced by a virtual photon with “field strength” 

equal to the matrix element of 


G a x 
let, +7,)v] - a! vy (1 svg} 


for creation of the neutrino pair. 

The author wishes to thank Professor William A. 
Fowler and Mr. M. Levine for suggesting a study 
of the reactions discussed here. 
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ABSTRACTS 


In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 








THE DIFFERENCE IN THE MULTIPLE SCAT- 
TERING OF ELECTRONS AND POSITRONS. 

B. P. Nigam, Department of Physics and Astron- 
omy, University of Rochester and The Research 
Division, Stromberg-Carlson, a Division of 
General Dynamics, Rochester, New York and 

V. S. Mathur, The Enrico Fermi Institute for 
Nuclear Studies, University of Chicago, Chicago, 
Illinois (Received November 9, 1960). 


The difference in the multiple scattering of 
electrons and positrons has been calculated on 
the basis of the work of Nigam, Sundaresan, 
and Wu. The results are compared with the ex- 
perimental work of Henderson and Scott and are 
found to be in good agreement. 


FREQUENCY FACTORS IN THE THERMALLY 
ACTIVATED PROCESS. R. Landauer and J. A. 
Swanson,* International Business Machines Re- 
search Center, Yorktown Heights, New York 
(Received November 3, 1960). 


Discussions of the rate, ve -U/kT for thermal- 
ly activated processes, are usually based on the 
phase-space distribution function for thermal 
equilibrium. Kramers has gone beyond this and 
for the particle in a bistable one-dimensional 
well has treated the relaxation to equilibrium as 
a Brownian motion problem in which the one- 
dimensional motion is coupled to a reservoir 
through a viscosity. Kramers’ arguments are 
readily extendable to many dimensions. In the 
overdamped case the reaction rate is reduced 
below the value derived from thermal equilibrium 
theory by the factor w,/n, where w, is the angu- 
lar frequency associated with the direction of 
steepest descent at the saddle point and 7 is the 
viscosity. In the underdamped case equilibrium 
theory is valid for many-dimensional systems, 
except for extreme degrees of underdamping. 


*Deceased. 
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DIELECTRIC RELAXATION DUE TO IMPURITY- 
VACANCY COMPLEXES IN NaCl CRYSTALS. 

R. W. Dreyfus, International Business Machines 
Research Center, Yorktown Heights, New York 
(Received September 19, 1960; revised manuscript 
received December 5, 1960). 


Sodium-ion vacancies are introduced into sodi- 
um chloride crystals by the addition of the fol- 
lowing divalent cation impurities: Mg*t+, Mn**, 
Ca++, Zn*+, Cd*+, andSr*++. Below 150°C, the 
virtual negative charge of a cation vacancy binds 
the majority of the vacancies to the divalent im- 
purity ions to form complexes. The reorientation 
of these dipolar complexes is observed in the 
present experiments as a transient dielectric 
polarization current occurring immediately after 
the application of a de electric field. Below 0°C, 
the current is made up of an initial exponential 
“fast polarization” due to complexes in which the 
vacancy is in the nearest sodium-ion site to the 
impurity ion, plus a “slow polarization” ascribed 
to vacancies jumping between more distant bound 
sites. The results for NaCl:Mn** crystals show 
that relaxation of vacancies which are nearest 
neighbors to the Mn** ion, occurs with an activa- 
tion energy of 0.680+0.010 ev. From the manner 
in which the relaxation time at a given tempera- 
ture varies with the radius of the impurity ion, 
the dominant jump mechanisms can be inferred. 
The major reason for the decrease in relaxation 
time with decreasing radius of the impurity ion 
appears to be the strong dependence of the jump 
rate of the impurity ion on ionic radius. 


EFFECTS OF APPLIED ELECTRIC FIELDS ON 
PARAMAGNETIC RESONANCE IN CHROME 
ALUMS. M. E. Browne, Lockheed Research 
Laboratories, Palo Alto, California (Received 
October 3, 1960; revised manuscript received 
December 7, 1960). 


The paramagnetic resonance spectra of single 
crystals of potassium chrome alum and ammo- 
nium chrome alum have been investigated in the 
presence of applied electric fields at 9.3 kMc/sec 
in the temperature range 4-300°K. The measure- 
ments were made in an effort to detect the pos- 
sible existence of a term of the form £y,0-E in 
the interaction Hamiltonian for an electron in an 
electromagnetic field. Such a term corresponds 
to the existence of an electron electric dipole 
moment and would result in a lifting of Kramers’ 
degeneracy for a magnetic ion in a crystalline 
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field. No such term was observed, and the 
upper limit on ~ obtained was é<0.007. A de- 
pendence of the dielectric loss on electric field 
was observed in both alums at microwave fre- 
quencies. 


ORDER- DISORDER MODEL THEORY FOR THE 
FERROELECTRIC EFFECT IN THE DIHYDRO- 
GEN PHOSPHATES. Michael E. Senko, Interna- 
tional Business Machines Research Center, 
Yorktown Heights, New York (Received Novem- 
ber 7, 1960). 


A term in polarization energy is added to the 
Slater-Takagi model of ferroelectricity in 
KH,PO,(KDP};). This free-energy function re- 
lates the shift in Curie point temperature due to 
deuteration of the compound with the measured 
change in saturation spontaneous polarization. 
The dielectric constant, €,, obtained by numer- 
ical differentiation of the free-energy function 
has a sharp maximum at the Curie temperature. 
With a change of parameters, the model also 
reproduces some of the details of antiferroelec- 
tric behavior. 


ISOTOPE EFFECT IN INTERMETALLIC DIF- 
FUSION. James G. Mullen,” Physics Depart- 

ment, University of Illinois, Urbana, Illinois 

(Received November 7, 1960). 


Ten measurements of the isotope effect for 
diffusion of Fe® and Fe™ in pure single crystals 
of silver and copper have been made over ap- 
proximately a 300-C° temperature interval be- 
low the melting point. The measured isotope 
effect was found to be about 2/3, significantly 
less than the value of unity expected from con- 
ventional reaction rate theory. 

It is shown from Vineyard’s extension of re- 
action rate theory that the measured isotope 
effect is to a good approximation a product of the 
Bardeen-Herring correlation factor and the frac- 
tion of the total translational kinetic energy asso- 
ciated with the diffusing tracer in a tracer- 
vacancy exchange. It is shown that the ring, 
interstitialcy, crowdion, and relaxation mechan- 
isms are not responsible for diffusion in the 
systems studied. 

The three-frequency theory of correlation for 
vacancy diffusion in fcc lattices is extended to 
show that the relative frequency factors are com- 


pletely specified by a knowledge of the correla- 
tion factor and the diffusion coefficients. An ex- 
pression for the difference in the Arrhenius Q 
and the activation energy due to the temperature 
dependence of the correlation factor is derived. 
This theory is shown to be inconsistent with the 
observed isotope effect if many-body effects are 
neglected. It is shown that a “long-range” re- 
pulsive interaction between the vacancy and im- 
purity would be necessary if the measured isotope 
effect is to be explained purely in terms of cor- 
relation. 


*Present address: Solid-State Science Division, 
Argonne National Laboratory, Argonne, Illinois. 


ELASTIC CONSTANTS OF AMMONIUM DIHY- 
DROGEN PHOSPHATE (ADP) AND THE LAVAL 
THEORY OF CRYSTAL ELASTICITY. Hans 
Jaffe, Clevite Corporation, Cleveland, Ohio and 
Charles S. Smith, Case Institute of Technology, 
Cleveland, Ohio (Received November 14, 1960). 


The acoustic shear wave velocities v 2 and 
vz, have been measured in ADP by the direct 
method of piezoelectric resonances, and by the 
direct pulse echo method. These velocities are 
found with each method to be equal within the 
experimental uncertainty. The resonance re- 
sults have a probable error of 0.1% for the dif- 
ference between the corresponding Laval elastic 
constants c,, andc,,. This experimental result 
disagrees with those previously obtained with 
ADP by indirect methods. 


PROTON MAGNETIC RESONANCE STUDY OF 
POLYCRYSTALLINE HCrO,. James A. Ibers, 
C. H. Holm, and C. R. Adams, Shell Develop- 
ment Company, Emeryville, California (Re- 
ceived October 21, 1960). 


A proton magnetic resonance study of poly- 
crystalline HCrO, as a function of magnetic field 
and temperature is presented. HCrO, is para- 
magnetic, and electron paramagnetic dipole as 
well as nuclear dipole effects lead to line broad- 
ening. The lines are asymmetric, and over the 
range of field 470 <H, < 9400 gauss and tempera- 
ture 77° <T < 294°K the asymmetry increases 
with increasing H, and decreasing T. An iso- 
tropic resonance shift of 0.017% to lower applied 
fields indicates a weak isotropic hyperfine con- 
tact interaction. The general theory of resonance 
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shifts is used to derive a general expression for 
the second moment M 2p» of a polycrystalline para- 
magnetic sample and is specialized to HCrO,. 

The theory predicts a linear dependence of Mo 

on [H,/(T +6), where @ is the experimentally 
determined Curie-Weiss constant. The experi- 
mental second moment M2p’ conforms to the re- 
lation 


Mop) = (3.96 + 0.06) + (0.084 + 0.002)[H,/(T +276)F, 


in agreement with theory. Hence the electron 
paramagnetic effects (slope) can be separated 
from the nuclear effects (intercept). The para- 
magnetic dipole effects provide some information 
on the particle shapes. The nuclear dipole effects 
provide some information on the motions of the 
hydrogen nuclei, but the symmetry of the O-H-O 
bond in HCrO, remains in doubt. 


OBSERVATIONS OF ELECTRON-HOLE CUR- 
RENT PINCHING IN INDIUM ANTIMONIDE. 
M. Glicksman and R. A. Powlus, RCA Labora- 
tories, Princeton, New Jersey (Received 
November 8, 1960). 


The temporal variation at high electric fields 
of the current and voltage in m -type InSb is ex- 
amined. The observations show the behavior 
characteristic of self-pinching of the electron- 
hole plasma current. There is good agreement 
with theory for the observed times for the pinch- 
ing process as a function of current. The cal- 
culated sum of the average energies of the elec- 
trons and holes in the plasma is 0.037 ev. Oscil- 
lations in the electric field of apparently differ - 
ent character are observed during the pinch, 
and when a longitudinal magnetic field is applied. 


EFFECTS OF A MAGNETIC FIELD ON HEAT 
CONDUCTION IN SOME FERRIMAGNETIC 
CRYSTALS. D. Douthett* and S. A. Friedberg, 
Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received November 10, 1960). 


The thermal conductivities of single-crystal 
bars of several ferrites have been measured 
over the temperature range 1.4° to 25°K both in 
zero external magnetic field and in a field of 
about 9.4 kilo-oersteds. For each specimen the 
zero-field conductivity is a smooth function of 
the temperature over the entire range considered. 
The magnitudes of the conductivities of the man- 
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ganese and cobalt-zinc ferrite specimens are 
quite low, approximately 1.5 mw/cm- °K at 4°K. 
Upon application of a magnetic field parallel to 
the direction of heat current flow, the conductivi- 
ties of the manganese ferrite and cobalt- zinc 
ferrite specimens are increased. The magneto- 
conductivities are rather large. For example, 
A4K/K,=0.2 for manganese ferrite at 2°K. Such 

a field produces no detectable change in the con- 
ductivity of the manganese-zinc ferrite. The ex- 
perimental results are interpreted on the basis 
of the assumption that the heat current is pri- 
marily a phonon current limited by phonon- 
magnon scattering. 


*Presently with Ramo-Wooldridge Research Labora- 
tory, Canoga Park, California. 


SUPERCONDUCTING TRANSITION IN ZINC AND 
GALLIUM. John F. Cochran* and D. E. Mapo- 
ther, Physics Department, University of Llinois, 
Urbana, Dlinois (Received November 14, 1960). 


The superconducting transitions of Ga and Zn 
have been studied to compare the supercooling 
of these elements with that shown by Al under 
similar conditions. Both Ga and Zn show con- 
siderably less tendency to supercool than Al 
despite their lower critical temperatures. The 
measurements also provide values of T, and 
(dH,./dT)T , for Ga and Zn for comparison with 
recently reported calorimetric data. Attention 
is called to the fact that discrepancies still exist 
between direct calorimetric determination of the 
specific heat discontinuity at the critical tem- 
perature and the value deduced thermodynamic- 
ally from critical-field data. It is shown that 
such discrepancies are conspicuously worse 
among the elements Zn, Ga, and Al than for the 
elements of higher T, (In, Sn, Ta, and Pb) 
where agreement is within experimental error. 


*Present address: Physics Department, Massachu- 
setts Institute of Technology, Cambridge, Massachu- 
setts. 


NUCLEAR QUADRUPOLE ABSORPTION IN IN- 
DIUM METAL. W. W. Simmons”* and C. P. 
Slichter, Department of Physics, University of 
Illinois, Urbana, Illinois (Received November 
7, 1960). 


The nuclear quadrupole resonance spectrum of 
indium metal has been observed and studied 
over the temperature range 4°K to 225°K. The 
four transition frequencies occur in ratios of 
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1:2:3:4, as expected from the crystal structure 
of indium. The lowest transition at 1.881 Mc/sec 
(at 4.2°K) gives a quadrupole coupling constant 

of (45.19+ 0.02) Mc/sec. 

The root second moment of 18 kc/sec is sub- 
stantially greater than the 4 kc/sec predicted 
from magnetic dipolar coupling among the nuclei, 
is temperature independent, but somewhat sam- 
ple dependent. It is shown that for a pure quad- 
rupole resonance, the pseudoexchange coupling 
produces an exchange broadening. This mecha- 
nism is shown to contribute about 10 kc/sec to 
the root second moment and predicts a Gaussian 
line shape, as observed. It is concluded that in 
some cases magnetic impurities must be pre- 
sent, but it is not known whether or not the nar- 
rowest line represents a natural width. 

The field gradient is computed assuming point 
ions and a uniform electron density by a tech- 
nique that converges rapidly with distance, but 
the answer is too small by a factor of 3. It is 
argued that it is important to consider the p va- 
lence electrons. The temperature dependence of 
frequency is very strong (a 23% decrease be- 
tween 4°K and 220°K), and unexplained. It appears 
to show that the field gradients cannot be com- 
puted assuming point ions and a smooth charge 
distribution. 

A resonance in the superconducting state is 
reported. 


*Now at Space Technology Laboratories, Inc., Los 
Angeles, California. 


PARAMAGNETIC RESONANCE HYPERFINE 
STRUCTURE OF TETRAVALENT Pa™ IN 
Cs,ZrCl,. J. D. Axe, Department of Chemistry 
and Lawrence Radiation Laboratory, University 
of California, Berkeley, California, and H. J. 
Stapleton and C. D. Jeffries, Physics Depart- 
ment, University of California, Berkeley, Cali- 
fornia (Received November 7, 1960). 


Paramagnetic resonance absorption at 3-cm 
wavelengths is observed for tetravalent Pa** 
in a single crystal of Cs,ZrCl, at helium tem- 
peratures. The observed spectra correspond to 
the allowed transitions (S,, 1, ~S,+1, I,) and 
the forbidden transitions (S,, I, +S,+1, I,* 1) 
of a system described by t the agin Samiitonton 
Hy = g8H- -S+Al- 8 -g, ‘BH: 1, with S=1/2, [=3/2, 
\A | /h =1578.6+1.4 Mc /sec, \g| =1.1423 + 0.0014, 
and |g,’|=8x10™. The errors indicate a small 
deviation from isotropy. It is further observed 


that g,,’/g<0, indicating that if &,' is positive, 

as is strongly indicated by the nuclear shell 
model, then g is negative. An additional electron- 
nuclear double resonance experiment is used to 
determine directly the nuclear magnetic moment 
u(Pa**")=1.96 nuclear magnetons. This value 
includes a correction of 9% due to perturbations 
of an excited state about 1900 cm™ above the 
ground-state doublet of Pa** in its octahedral 
crystal field. A lower frequency double-reso- 
nance experiment is used to measure the weak 
hyperfine interaction of the Pa‘* ion with its Cs'*? 
neighbors, of order A’’/h=0.5 Mc/sec. 


INTERACTION OF SLOW ELECTRONS WITH IN- 
SULATING CRYSTALS. I. COMPARISON OF 
ELECTRON AND PHOTON ABSORPTION COEF- 
FICIENTS FOR KCl AND KBr. William J. Fre- 
dericks and Charles J. Cook, Stanford Research 
Institute, Menlo Park, California (Received 
October 28, 1960). 


The electron absorption coefficient for pure 
KCl and additively colored KCl and KBr was 
measured in the energy range from 0.2 to 12 ev. 
These data are compared with photon absorption 
measured in the same crystals in the energy 
range from 1.2 to 6.5 ev and with published opti- 
cal and photoemission data above 6.5 ev. Cer- 
tain imperfection interactions were found at 
corresponding photon and electron energies. 

The surface transition layer, even for pure KCl, 
was observed to contain F and V, centers. 

There was no evidence of direct electron- exciton 
interactions. The valence band of KCl was ob- 
served to be 8.6+0.2 ev below the vacuum level; 
this agrees with Eby, Teegarden, and Dutton’s 
measurements. Upon electron bombardment, an 
optical absorption near the first fundamental 
band was found to develop. Its characteristics 
do not coincide with any of the absorption bands 
previously reported in this energy region. 


VISIBLE LUMINESCENCE OF RARE EARTH- 
YTTRIUM GALLIUM GARNETS. S. P. Keller 
and G. D. Pettit, International Business Ma- 
chines Research Center, Yorktown Heights, New 
York (Received October 3, 1960). 


Yttrium gallium garnet (Y,Ga,O,,) has been 
prepared with small percentages of different 
rare earths substituted for yttrium. Garnets 
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activated with Pr, Sm, Eu, Tb, Dy, Ho, Er, 
and Tm are all luminescent. The emission and 
excitation spectra of the samples have been 
measured at 77°K. The data have been analyzed 
in terms of the atomic energy levels of the im- 
purity ions and the effects of crystalline field 
and of phonon interactions. Wherever possible, 
comments are made about whether the crystal 
field can be treated as possessing cubic or low- 
er symmetry. 


ELECTRICAL RESISTIVITY OF LANTHANUM, 
PRASEODYMIUM, NEODYMIUM, AND SAMAR- 
IUM. J. K. Alstad, R. V. Colvin,* S. Legvold, 
and F. H. Spedding, Institute for Atomic Re- 
search and Department of Physics, Iowa State 
University, Ames, Iowa (Received November 9, 
1960). 


The electrical resistivities of polycrystalline 
samples of La, Pr, Nd, and Sm are reported 
in the temperature range 1.3°K to 300°K. La 
exhibits a superconducting transition at 5.8°K. 
The curve for Pr has slope changes at 61°K and 
95°K. The Nd curve shows small jumps at 5°K 
and 20°K. Sm shows slope changes at 14°K and 
106°K. 

*Now at the Edgar C. Bain Laboratory of the U. S. 
Steel Corporation, Monroeville, Pennsylvania. 


LOW-TEMPERATURE MAGNETIC SUSCEPTI- 

BILITIES OF THE HYDRATED ACETATES OF 

Ni*+, Co**, AND Mn**. R. B. Flippen* and 

S. A. Friedberg, Carnegie Institute of Technol- 
ogy, Pittsburgh, Pennsylvania (Received August 
9, 1960). 


The magnetic susceptibilities of powdered 
Ni(C,H,O,),°4H,O, Co(C,H,O,),-4H,O, and 
Mn(C,H,O,),°4H,O have been measured in the li- 
quid hydrogen and helium ranges by an ac mutual 
inductance method. Measurements on the cobalt 
salt were extended to 0.4°K, revealing a sharp 
maximum at 0.6°K. A similar, more pronounced 
susceptibility peak is observed in the manganous 
acetate at 3.2°K. Single-crystal measurements 
have been made on the latter salt along the mono- 
clinic b axis as well as in the ac plane in external 
fields up to 1000 gauss. Saturation and hysteresis 
are seen in the a direction associated with the ap- 
pearance of a weak ferromagnetic moment below 
3.18°K. Mn(C,H,O,),-4H,O appears to be a weak 
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ferromagnet of the kind discussed recently by 
Dzyaloshinsky and Moriya. 


*Now at Central Research Department, Experi- 
mental Station, E. I. duPont de Nemours and Co., 
Wilmington, Delaware. 


BEHAVIOR OF EXCITED ELECTRONS AND 
HOLES IN ZINC SULFIDE PHOSPHORS. Shigeo 
Shionoya,* Hartmut P. Kallmann, and Bernard 
Kramer,t Department of Physics, New York 
University, Washington Square, New York, New 
York (Received October 11, 1960). 


The luminescent properties of silver- and 
copper -activated ZnS phosphors have been investi- 
gated at room and liquid air temperatures using 
ultraviolet, blue light, and 8-ray excitation. In 
many cases it was found that the deficiency area 
was considerably smaller than the glow-curve 
area, especially for 8-ray excitation of ZnS:Cu. 
This is explained by the hole production causing 
quenching during §-ray excitation, while during 
the glow curve the smaller amount of holes pro- 
duced causes a much smaller amount of quenching. 
However, even with blue light where excitation 
does not directly produce holes, the deficiency 
area is smaller than the glow-curve area. 

The effect of infrared radiation on the fluo- 
rescence, phosphorescence, and glow curve has 
been investigated, and the quenching and stimu- 
lation at room and liquid air temperature re- 
ported. For ZnS:Cu the long-wavelength infrared 
band (~12 000 A) produces quenching at room 
temperature but only stimulation at low tempera- 
ture. 

A determination of trap density from the glow 
curves shows ~10'’/cm® for ZnS:Cu and ~ 10'®/cm* 
for ZnS:Ag at liquid air temperatures. 

Finally, a theoretical discussion of the possible 
transition mechanisms giving rise to these effects 
is given. 


*Work performed while on leave from Yokohama 
National University, Yokohama, Japan. Presently 
associated with The Institute for Solid State Physics, 
University of Tokyo, Tokyo, Japan. 

TtAlso Department of Physics, Hunter College, New 
York, New York. 


GREEN’S FUNCTION METHOD FOR QUANTUM 
CORRECTIONS TO THE THOMAS- FERMI MO- 
DEL OF THE ATOM. Gene A. Baraff and Sidney 
Borowitz, Physics Department, New York Uni- 
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versity, New York, New York (Received August 
12, 1960; revised manuscript received Decem- 
ber 13, 1960). 


A systematic method is presented for deriving 
the Thomas- Fermi equation for an atom and the 
quantum corrections from the many-body de- 
scription. The novel feature of the method is 
that it does not require any a priori assumptions 
about the assignment of electrons to fully occu- 
pied single-particle states or about the distribu- 
tion of electrons in phase space, but shows in- 
stead that the distribution which is usually as- 
sumed, or derived from the assumption of fully 
occupied single-particle states, is a direct con- 
sequence of specifying that the many- particle 
system is in its ground state. The procedure 
used in the derivation is the expansion of the 
mixed position- momentum representation of the 
Green’s function in a series of powers off. The 
lowest order term is found to correspond with 
the Thomas- Fermi density. The form of the 
higher order terms, which are to be considered 
as corrections to zeroth order term, depends on 
the approximations made in the many-body equa- 
tions for obtaining the Green’s function. This 
paper deals only with the Hartree- Fock approxi- 
mation, but the methods presented here allow 
generalization to other approximations which can 
include correlation effects. 


ENERGY DISTRIBUTION OF ELECTRONS FROM 
IONIZING COLLISIONS OF HEAVY PARTICLES. 
H. W. Berry, Syracuse University, Syracuse, 
New York (Received October 17, 1960). 


The energy distribution of electrons produced 
in ionizing collisions of ions and fast atoms with 
atoms of the parent gas at rest has been meas- 
ured for A and He. The energy of the incident 
particles ranged from 0.30 to 3.0 kev. Electrons 
released at 90° with the incident ion or atom 
beam were analyzed in an electrostatic energy 
selector which consisted of 90° segments of 
coaxial cylinders. The use of a fine wire grid in 
the region in which the ionization occurred allow- 
ed a positive differentiation to be made between 
the ionization electrons and the secondary elec- 
trons from any metal surfaces. The electron 
distributions for argon ions and atoms in argon 
are similar and show a rapid decrease in the 
yield as the energy increases from zero to about 
4 ev followed by a plateau and a flat maximum 
near 12 ev. The distributions for He also show 





the initial rapid decrease continuing to a small 
maximum which for the ion is at about 31 ev and 
for the neutral atom at about 16 ev. These dis- 
tributions are almost collision energy independ- 
ent, though with an increase in the over-all yield 
with increasing energy. 


ISOTOPE SHIFTS IN PALLADIUM. R. H. Hughes 
and F. A. Sharpton,* Physics Department, Uni- 
versity of Arkansas, Fayetteville, Arkansas 
(Received November 14, 1960). 


Isotope shifts in the 4d° 5s 'D,—4d* 5p °F,° 
transition at \4212 A in the first spectrum of Pd 
have been studied with the use of enriched iso- 
topes. The shifts are similar to those found by 
Kuhn and Warner in the 4d° 5s °D,—4d® 5p °F ,° 
transition at \3405 A. The most interesting fea- 
ture is the maximum in the isotope shift which 
appears at the neutron number pair 60-58. 


*Now at Hendrix College, Conway, Arkansas. 


BOUNDS ON SCATTERING PHASE SHIFTS FOR 
COMPOUND SYSTEMS. Leonard Rosenberg and 
Larry Spruch, Physics Department, Washington 
Square College and Institute of Mathematical 
Sciences, New York University, New York, New 
York (Received November 9, 1960). 


An extension of recently developed methods 
determines a rigorous upper bound on (-k coty)~*, 
where 7 is the phase shift, for the general one- 
channel scattering process. The method unfor- 
tunately requires truncation of the various po- 
tentials, but it should generally be possible, in 
practice, to so truncate the potentials that the 
difference between the phase shifts of the origi- 
nal problem and of the problem for which a 
bound is obtained is insignificant. 

In the course of the development it is necessary 
to introduce, for compound-system scattering, 
an absolute definition of the phase shift, not 
simply a definition modulo 7. The definition 
chosen is to take the projection of the full scat- 
tering wave function on the ground-state wave 
function of the scattering system, and to treat 
the resultant one-coordinate wave function as if 
it were the scattering wave function for a par- 
ticle on a center of force. Though irrelevant 
with regard to the determination of a bound on 
cot7, it is interesting that at least for some 
simple cases this definition automatically in- 
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creases the phase shift by at least 1 whenever 
the Pauli principle introduces a spatial node into 
the scattering wave function. The triplet scatter- 
ing of electrons by hydrogen atoms provides an 
example. 


ISOTOPE SHIFT IN THE ARC SPECTRUM OF 
NICKEL. Daniel J. Schroeder and J. E. Mack,* 
University of Wisconsin, Madison, Wisconsin 
(Received September 19, 1960). 


The isotope shift in 31 spectral lines in the 
nickel arc spectrum has been determined by the 
use of a Fabry-Perot interferometer. The nor- 
mal mass shifts were calculated (=+ 0.025 cm™! 
between Ni® and Ni™) and subtracted from the 
observed isotope shifts. The differences were 
attributed to the specific mass and field effects. 
The relative shifts of levels of four configura- 
tions were deduced from the observed line shifts, 
these being the “complex” configurations, 3d*4s? 
and 3d*4s4p, and the two-electron configurations, 
3d°4s and 3d°4s. It was shown that the shifts due 
to the specific mass effect are a significant part 
of the observed shifts. Perturbations due to in- 
terconfiguration interactions were postulated to 
explain some of the observed shifts. The isotope 
shift to be expected between Ni®® and Ni™ on the 
basis of field-effect calculations is about -0.020 
cm™' for a single 4s electron, while the shifts 
observed are as large as +0.190 cm™'. A large 
fraction of this shift must therefore be attributed 
to the specific mass effect. By noting the devia- 
tions of the relative shifts between adjacent pairs 
of even isotopes from those predicted by mass- 
effect theory, it was possible to deduce the re- 
lative field effect. The relative level shift re- 
sulting from the field effect is nearly the same 
for the adjacent isotope pairs, 60-62 and 62-64, 
while the relative level shift for the isotope pair, 
58-60, is approximately 0.004 cm™ larger than 
that for the other adjacent isotope pairs. The 
arrangement of neutrons in the outermost nuclear 
shells is believed to account for this difference. 
Within the experimental error the level shift of 
the Ni®™ relative to the neighboring even isotopes 
is such that there is no odd-even staggering of 
the levels. 


*Currently on leave as Scientific Attaché, Stockholm, 
Sweden. 
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ANALYSIS OF GAMMA-GAMMA POLARIZATION- 
DIRECTIONAL CORRELATIONS INVOLVING 
MULTIPOLE MIXTURES. M. B. Martin, R. G. 
Arns, and M. L. Wiedenbeck, Harrison M. Rand- 
all Laboratory of Physics, The University of 
Michigan, Ann Arbor, Michigan (Received Au- 
gust 19, 1960). 


A method is described for the graphical anal- 
ysis of the gamma-gamma polarization-direc- 
tional correlation when one or both gamma rays 
involve a mixture of dipole and quadrupole radia- 
tion. The case in which the polarization of the 
mixed gamma ray is observed is treated in de- 
tail. 


STATISTICAL THEORY OF GAMMA-RAY 
SPECTRA FOLLOWING NUCLEAR REACTIONS. 
E. S. Troubetzkoy, Nuclear Development Cor- 
poration of America, White Plains, New York 
(Received October 18, 1960). 


A theory predicting y-ray spectra following 
intermediate or high-energy nuclear reactions 
is derived on the basis of the statistical theory. 
The equations take a particularly simple form 
if the assumption is made that all radiative 
transitions are of the electric dipole type. The 
theory is applied to two specific reactions: in- 
elastic neutron scattering and thermal neutron 
capture. Numerical calculations of spectra 
arising from thermal neutron capture by two 
gadolinium isotopes are shown to compare well 
with experiments. 


NEUTRON-CAPTURE GAMMA-RAY SPECTRA 
OF THE NICKEL ISOTOPES. P. A. Treado and 
P. R. Chagnon, Physics Department, The Uni- 
versity of Michigan, Ann Arbor, Michigan (Re- 
ceived September 14, 1960). 


The spectra of gamma radiation following pile- 
neutron capture in separated isotopes of nickel 
have been measured with a three-crystal scintil- 
lation pair spectrometer as well as with a single- 
crystal scintillation counter. Energies and in- 
tensities of the gamma rays from each isotope 
are reported. Each of the lines observed with 
natural nickel has been assigned to Ni™(n, y)Ni*®, 
Ni™(n, y)Ni™, or Ni®(n, y)Ni®. In some cases, 
gamma rays of nearly identical energies are 
found in two, or in all three, isotopes. Additional 
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gamma rays not resolved in the natural nickel 
spectrum appear in Ni®(n, y)Ni®® at 7.65 Mev, 

and in Ni®(n, y)Ni™ at 4.70 and 5.55 Mev. The 
spectrum of Ni®(n, y)Ni® is dominated by a strong 
line at 6.80 Mev, presumably the ground-state 
transition. 


ELASTIC SCATTERING OF ALPHA PARTICLES 
BY N’®. H. Smotrich, K. W. Jones,* L. C. 
McDermott, and R. E. Benenson,t Columbia 
University, New York, New York (Received 
August 24, 1960). 


Absolute differential cross sections for the 
elastic scattering of alpha particles by N*® have 
been measured in a differentially-pumped gas 
scattering chamber. The measurements were 
made at center-of-mass angles of 169.1, 149.5, 
140.8, 125.3, 90.0, and 70.0 degrees for alpha- 
particle energies from 1.75 to 5.50 Mev, cor- 
responding to 5.37- to 8.33-Mev excitation of 
the compound nucleus F**. The experimental 
widths of the levels observed below approxi- 
mately 3.5-Mev bombarding energy are general- 
ly narrow, in most cases less than 10 kev. Above 
3.5 Mev a marked increase in the level width was 
observed. As a result of an analysis based on 
the Wigner-Eisenbud reaction theory, values of 
J, 7, Ey, and y,* were assigned to 16 levels in 
yp. 

*Now at the Ohio State University, Columbus, Ohio. 

tPermanent address: City College of New York, New 
York, New York. 


ANALYSIS OF LONG-RANGE ALPHA- EMISSION 
DATA. R. D. Griffioen and J. O. Rasmussen, 
Lawrence Radiation Laboratory, University of 
California, Berkeley, California (Received No- 
vember 2, 1960). 


Alpha barrier penetrabilities for the long-range 
alphas of Po?!” and Po?“ are calculated by using 
the diffuse exponential nuclear potential derived 
from optical- model analysis of alpha elastic- 
scattering data. The calculations are made on 
the same basis as reported by Rasmussen in two 
previous publications. Partial half-lives for 
alpha and gamma emission are calculated on the 
assumption that the long-range alpha decay is 
unhindered with respect to the ground-state alpha 
decay. 


METHOD TO OBTAIN THE ENERGY DEPEND- 
ENCE OF REACTION CROSS SECTIONS WHEN 

A COMPOUND NUCLEUS IS FORMED. Raymond 
Fox and Richard D. Albert, Lawrence Radiation 
Laboratory, University of California, Livermore, 
California (Received November 10, 1960). 


It is found that measuring the energy distri- 
butions of emitted particles at two different in- 
cident energies and two reaction cross sections 
is sufficient to obtain the energy dependence of 
a reaction cross section when it is known that a 
compound nucleus is formed. The cross sections 
obtained are a function of the excitation energy 
of the residual nucleus. They are the usual re- 
action cross sections at zero excitation energy 
when the assumption is made that the reaction 
cross section does not vary with the excitation 
energy of the target nucleus. By comparing the 
reaction cross sections thus obtained with meas- 
ured nonelastic cross sections, the compound 
elastic cross section may also be obtained. An 
example using this method is given in which two 
measured energy distributions and two measured 
cross sections were used to obtain the energy de- 
pendence of the reaction cross section of nickel 
in the 2.5- to 7-Mev energy range. 


ENERGY LEVELS IN Zr™ EXCITED BY THE 
(d, p) REACTION. R. L. Preston,” H. J. Mar- 
tin, Jr.,! and M. B. Sampson, Indiana University, 
Bloomington, Indiana (Received November 3, 
1960). 


Eleven energy levels in Zr™ were excited. 
The excitation energies of the levels are 0, 1.22, 
2.07, 2.19, 2.56, 2.84, 3.05, 3.25, 3.45, 3.65, 
and 3.87 Mev. Angular distributions were meas- 
ured for the first four levels and assignments 
were made using the shell model and the Butler 
stripping theory. The assignments are 2d,, for 
the ground state, 3s,,. for the 1.22-Mev state, 
2d,,. for the 2.07-Mev state, and 1g,, for the 
2.19-Mev state. The application of the Butler 
theory to the data gave reasonable and useful 
fits. The ground-state @ value was measured 
and is 5.02+0.03 Mev. 


*Present address: 35 Bedford Park Avenue, Rich- 
mond Hill, Ontario, Canada. 

TNow on leave at the Brookhaven National Laboratory, 
Upton, New York. 
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NATURAL ALPHA RADIOACTIVITY IN MEDIUM- 
HEAVY ELEMENTS. Ronald D. Macfarlane* 

and Truman P. Kohman, Department of Chemis- 
try, Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received November 14, 1960). 


A large cylindrical ionization counter accomo- 
dating samples up to 1200 cm? in area has been 





Alpha particle 


Specific activity of 


used for measurements of natural alpha radio- 
activity in medium-heavy elements. Low-back- 
ground techniques and multichannel pulse analy- 
sis are employed. The method has greater energy 
resolution and yields better counting statistics 
than the nuclear emulsion technique, but does 

not have as great a sensitivity. The results ob- 
tained from measurements on natural elements 
and isotopically enriched samples are: 


Half-life or limit 





energy natural element for nuclide 
Nuclide (Mev) (dis sec g-') (years) 
Cel#2 vee <0. 0002 >5 x10 
Na!44 1.83 +0. 02 0.0093 +0.0011 (2.4 +0.3) x10 
Sm'46 ee <0.013 <3 x10° 
Sm‘4? 2.23 +0. 02 116 +5 (1. 15 +0. 05) x10"! 
Sm'“8 tee <0. 04 >2 x10"4 
Sm'9 see <0.01 >1 x10 
Ga's? 2.14 +0. 03 0. 00156 +0. 0010 (1. 08 +0. 08) x10"* 
Hf!™ 2.50 +0. 03 0. 000065 +0. 000013 (2.0 +0. 4) x10% 
wise see <0. 00013 >9 x10!4 
pr? 3.11+0.03 0.012 +0. 001 (6.9 +0. 5) x10!! 
Hg! tee <0. 0007 >1 x10"4 


*Present address: Lawrence Radiation Laboratory, University of California, Berkeley, California. 


PREDICTIONS OF SPONTANEOUS FISSION 
HALF-LIVES FOR HEAVY NUCLEI. David W. 
Dorn, Lawrence Radiation Laboratory, Univer- 
sity of California, Livermore, California (Re- 
ceived November 4, 1960). 


Swiatecki’s work on correlation of spontaneous 
fission half-lives has been modified and extended 
to include elements beyond Z=100. The values 
of spontaneous fission half-lives predicted on 
these bases are unexpectedly high. For example, 
the partial half-life for Z =106, A =271 is pre- 
dicted to be about 13 years. 


DECAY OF 49-min Cd"* AND 5.1-sec In**. 
Chester E. Gleit* and Charles D. Coryell, De- 
partment of Chemistry and Laboratory for Nu- 
clear Science, Massachusetts Institute of Tech- 
nology, Cambridge, Massachusetts (Received 
November 2, 1960). 


The decay properties of 49-min Cd""® and its 
daughter 5.1-sec In™* are established. The 
Cd'"* has a low 8-decay energy and seems to be 
an allowed transition. The In'!® has (4.2 + 0.4)- 
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Mev £-decay energy, with about 15% of the de- 
cays going to the 1.22- Mev excited state of Sn™"*, 
both £ transitions being allowed. It is proposed 
that a previously identified 4.5- min In species 
is a high-spin isomer In"*” not formed in the 

B decay of Cd""* or formed directly in uranium 
fission. 


*Present address: Westinghouse Electric Corpora- 
tion, Bettis Atomic Power Laboratory, Pittsburgh, 
Pennsylvania. 


DECAY OF Ne’®. J. W. Butler and K. L. Dunning, 
Nucleonics Division, U. S. Naval Research La- 
boratory, Washington, D. C. (Received Novem- 
ber 2, 1960). 


The radioactive nuclide Ne** has been produced 
in the O'*(He*, n)Ne*® reaction by 5.2-Mev He® 
particles from the NRL 5- Mv Van de Graaff 
accelerator. The target, CaO on a Pt backing, 
was about 300 kev thick to the incident beam. 
The gamma rays following the decay of Ne*® 
were detected by a 3-in. diam by 3-in. Nal(T1) 
crystal and a 256-channel pulse-height analyzer. 











VOLUME 6, NUMBER 2 


PHYSICAL REVIEW LETTERS 


January 15, 1961 





Only one nuclear gamma ray was observed. Its 
energy, measured with respect to the Na”* gam- 
ma ray at 1.27364 0.0016 Mev, is 1.041+ 0.005 
Mev. Therefore, this gamma ray is assigned to 
the 1.04-Mev state of F’®. This state is con- 
sidered to be the isobaric spin analog state of 
the Ne’® ground state, and therefore has quantum 
parameters 0*, T=1. The branching to this ex- 
cited state was measured to be 0.07+ 0.02 (in- 
cluding a calculated contribution of 0.00006 from 
orbital electron capture), and that to the ground 
state of F’*® was measured to be 0.934 0.02 (in- 
cluding a calculated contribution of 0.0002 from 
orbital electron capture). The ft values are, 
respectively, 3030+ 880 and 1120+ 70 sec. The 
half-life of Ne*® was measured to be 1.46+ 0.07 
sec. 


O** DECAY ENERGY AND THE FERMI INTER- 
ACTION CONSTANT. J. W. Butler and R. O. 
Bondelid, Nucleonics Division, U. S. Naval Re- 
search Laboratory, Washington, D. C. (Received 
November 10, 1960). 


The threshold energy of the C’*(He*,)O” re- 
action has been measured precisely with the use 
of a Van de Graaff accelerator and an electro- 
static analyzer, the value being 1436.2+0.9 kev. 
From this value, the O'* beta- decay end-point 
energy (for decay leading to the 2.312- Mev state 
in N**) is computed to be 1.8000+ 0.0065 Mev, 
based on the 1956 table of masses, and 1.8097 
+ 0.0015 Mev, based on the 1960 table of masses. 
A revised value of the Fermi interaction constant 
in beta decay is calculated and applied in the 
conserved-vector-current theory of Feynman and 
Gell-Mann. When (1) the radiative corrections 
and other corrections are applied to the decay of 
Oo“, (2) the corrected ft value is used to com- 
pute the vector coupling constant in beta decay, 
(3) the value of this vector coupling constant is 
assumed to be the same as that for the muon de- 
cay and is used to calculate the lifetime of the 
muon, and (4) this lifetime is corrected for ra- 
diation effects, the predicted mean life of the 
muon becomes 2.289+ 0.013 psec (based on the 
1960 table of masses and the radiative correc- 
tions of Kinoshita and Sirlin) or 2.245 0.013 
usec (based on the 1960 table of masses and the 
corrections, radiative and otherwise, of Durand 
and collaborators). The former value is 3.6% 
greater than a weighted average of several re- 
cent measurements of the muon mean life, 2.210 


+0.003 sec, while the latter value is only 1.6% 
greater, and is within the combined experimen- 
tal and theoretical uncertainties. However, the 
definitions of coupling constants used by Durand 
and collaborators differ somewhat from those 
used by others. 


PHOTOPRODUCTION AND DETECTION OF THE 
TWO-MESON BOUND STATE. Jack L. Uretsky 
and T. R. Palfrey, Jr., Purdue University, La- 
fayette, Indiana (Received August 3, 1960; re- 
vised manuscript received December 14, 1960). 


We investigate the possibility of the detection 
of a photoproduced two-pion bound system. The 
general detection problem is discussed briefly; 
the branching ratio between 27° and 2y° decay 
modes is calculated; and the total photoproduc- 
tion cross section is estimated in terms of the 
binding energy of the two-pion state both by field- 
theoretic and by phase-space arguments. We 
conclude that if the binding energy is of the order 
of 10 Mev the state should be detectable in photo- 
production experiments, and the binding energy 
should be measurable. 


EXCITATION FUNCTIONS OF (p, 2 NUCLEON) 
REACTIONS. P. P. Strohal* and A. A. Caretto, 
Jr., Department of Chemistry, Carnegie Insti- 
tute of Technology, Pittsburgh, Pennsylvania 
(Received September 26, 1960; revised manu- 
script received December 5, 1960). 


Excitation functions for the (p,pn) reaction 
have been determined in the energy range 250 
Mev to 440 Mev for the target nuclei As™, Br®, 
zZr®, ['?7, Ce*#?, w'®*, and Re’®’. Both isomer 
products were determined in the case of the 
Br®(p,pn)Br® reaction. Also, excitation func- 
tions were determined for the (p, 2p) reaction on 
Zr®, Ce’, and W’™ and for the (p, 2m) reaction 
on Zr®™, All cross sections were measured re- 
lative to the Al*?’(p, 3pn)Na™ excitation function. 
The (p,pn) excitation functions are only moder- 
ately sensitive to energy, showing an increase by 
about 30% for iodine, 20% for arsenic and cer- 
ium, and a 12% decrease for zirconium. The 
remaining excitation functions are flat. The 
ratio of the OF, to Op 2p is 1.62+0.09, for the 
three targets measured at 440 Mev. The various 
mechanisms for the production of reactions of 
this type are examined in the light of present 
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understanding. The suggestion is made that, 
despite the relative insensitiveness of these ex- 
citation functions to energy, (p,pn) and (), 2p) 
reactions may take place by more than one 
mechanism in this energy region. A comparison 
is made with the results of recent Monte Carlo 
nuclear cascade calculations applied to calcula- 
tion of (p, pn) reaction excitation functions. 


*Present address: Institut Ruder Boskovic, Zagreb, 
Yugoslavia. 


RADIATIVE DECAY MODES OF THE MUON. 
R. R. Crittenden and W. D. Walker, University 
of Wisconsin, Madison, Wisconsin, and J. Bal- 
lam, Michigan State University, East Lansing, 
Michigan (Received August 29, 1960). 


A 5-in. Freon bubble chamber was used to 
search for the following decay modes of the ».* 
meson: 


(1) ut ~et +Y; 

(2) ut++et+e~+e*, 

(3) pt +et+vy° +P +y, 

(4) utwet+y° +P +et+e~. 
Two exposures were made at the Carnegie Insti- 
tute of Technology synchrocyclotron. A total of 
200 000 pictures were taken, yielding 3.3 x10° y*- 
meson decays. 

A total of 3x10° u* decays were examined for 
mode (1). No decays consistent with this mode 
were found. The upper limit on the branching 
ratio Rygq was found to be 

R_ ,=(u* wet +y) Apt wet +9 +0) <2.5x1075. 
rad 

A total of 3.3x10° u* decays were scanned for 


mode (2) and no such decays were observed. The 
limit on the branching ratio Rg, was found to be 


Re, =(u +e +e +e °)/(p* +e" +°+0)<4x10™. 


The internal bremsstrahlung rate (mode 3) was 
measured for two values of Eg (the minimum 
photon energy detected). The results were 


Ryp=(u* et ++ +y) Aut er +y° +P); 


= ~2 = . 
Rip (1.4+0.4)x107?, BE, 10 Mev; 


0 


= “3 = 
Ryp=(3.341.3)x10", E, .=20 Mev. 
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The rate of internal conversion of internal 
bremsstrahlung (mode 4) was found to be 


Rig = (hi wet ++ +e* +e )/(u* e* + +P) 


=(2.2+1.5)x1075, £,=10 Mev, 


where E,, is the minimum energy of the internally 
converted y ray. 

A summary is given of previous experiments 
on these decay modes and results are discussed 
with special reference to the intermediate- 
boson scheme of weak four-fermion interactions. 


CLOUD-CHAMBER STUDY OF HARD COLLI- 
SIONS OF COSMIC-RAY MUONS WITH ELEC- 
TRONS. R. F. Deery* and S. H. Neddermeyer, 
Department of Physics, University of Washing- 
ton, Seattle, Washington (Received November 1, 
1960). 


The energy distribution of the secondary elec- 
trons produced in targets of carbon or paraffin 
by the yp-e scattering process for muon momenta 
in the range 5-50 Bev has been measured for 
electron energies up to 20 Bev, or c.m. momen- 
tum transfer up to 140 Mev. A vertical array of 
three cloud chambers immersed in a magnetic 
field of 11000 gauss was used with a fourfold 
coincidence system. Two flat rectangular pro- 
portional counters, suitably biased, together 
with two Geiger- Muller trays, provided fair re- 
jection of uneventful penetrating particles and a 
high efficiency for selection of the narrow elec- 
tronic showers characteristic of the high-energy 
electromagnetic events. On 5900 counter-trig- 
gered photographs there were 291 accepted 
events, having one or more (+) electrons with 
energy 20.10 Bev, believed to originate in p-e 
collisions in the (carbon or paraffin) target above 
the top chamber, from incident muons in the 
momentum interval 5-50 Bev. The data are com- 
pared with a calculation based on the Bhabha for- 
mula for spin 1/2 muons, taking into account the 
momentum distribution of the incident muons, 
the energy loss and shower development in the 
target and the chamber walls, and a theoretical 
efficiency factor. Arguments are given to show 
that direct pair production and bremsstrahlung 
of the muons in the target and in the Pb shield 
above the apparatus produce negligible effects. 
The experiment permits a reliable measurement 
of only the relative distribution. When arbitrarily 
normalized, the calculated distribution is in 
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fairly good agreement with the data, except for 

a small systematic difference suggesting an ex- 
cess of observed events for the harder collisions. 
Although the discrepancy is interpretable as a 
statistical fluctuation, the data are fitted much 
better over the entire range when the basic 
cross section is modified by a “form factor,” 

F? (greater than unity), with F =1+1q@lA,?, 
where gq is the invariant of the 4-momentum 
transfer in units of k, and A, is the Compton 
wavelength of the muon. This may be the first 
indication of a deviation from standard quantum 
electrodynamics for hard p-e collisions. More 
strongly it shows that, if there is a deviation, 

it is not representable by a form factor less than 
unity. 


*Now at the California Institute of Technology, 
Pasadena, California. 


SPIN IN CLASSICAL AND QUANTUM THEORY. 
H. C. Corben, Space Technology Laboratories, 
Los Angeles, California (Received March 31, 
1960; revised manuscript received October 13, 
1960). 


The classical equations of motion of a charged 
point particle with intrinsic spin under the in- 
fluence of an external electromagnetic field are 
restated and compared with the Heisenberg equa- 
tions of motion derived from the Dirac theory. 
The partition of angular momentum between par- 
ticle and field in the classical theory is con- 
trasted to the Dirac theory of electron spin. The 
analogy between the Dirac equation and the theory 
of parametric amplification is pointed out. 

A free spinning point particle moving accord- 
ing to the laws of classical relativistic point- 
particle mechanics may move along a helix. 

The sum of the intrinsic spin o and the angular 
momentum of the helical motion in this classical 
analog of Zitterbewegung is an effective spin 
vector S which is a constant of the motion. Be- 
cause of this internal motion, the effective mass 
M of the particle differs from the mass m which 
is ascribed to it in the equations of motions. 
Solutions are found in which S is parallel or 
antiparallel to the momentum, and the sign of M 
is determined by the helicity. When placed ina 
uniform electromagnetic field, the particle be- 
haves as if it had a rest mass M and a magnetic 
moment e0/Mc, in addition to any explicit mag- 
netic moment that may be ascribed to it. 


CONTINUATION OF SCATTERING AMPLITUDES 
AND FORM FACTORS THROUGH TWO- PARTI- 
CLE BRANCH LINES. Reinhard Oehme, Enrico 
Fermi Institute for Nuclear Studies and Depart- 
ment of Physics, University of Chicago, Chicago, 
Illinois (Received November 3, 1960). 


It is shown that scattering amplitudes and form 
factors have two-particle branch lines which 
connect two Riemann sheets. For partial wave 
amplitudes and form factors the dispersive parts 
and, except for square root factors, the absorp- 
tive parts are regular functions in the cut energy 
plane except for isolated poles, physical inelas- 
tic cuts, and left-hand branch lines. In order to 
show this it is assumed that, for particles with- 
out composite structure, the amplitudes have 
only such singularities in the physical sheet 
which correspond to absorptive processes. The 
analytic properties of absorptive parts are used 
for a general discussion of structure singulari- 
ties (anomalous thresholds). It is shown that 
these structure cuts are extensions of left-hand 
branch lines. An example is given of a disper- 
sion relation on the Riemann surface in which 
the integral over the two-particle branch line is 
eliminated. 


RELATIVISTIC PARTICLE DYNAMICS. B. Ba- 
kamjian, The Pennsylvania State University, 
University Park, Pennsylvania (Received Novem- 
ber 10, 1960). 


It is shown that an instant form of relativistic 
particle dynamics can be constructed which dis- 
plays the symmetry properties inherently pre- 
sent in the point form. In this new instant form, 
interaction terms are introduced in the energy- 
momentum four-vector of the system; physically 
this is more justifiable than the more customary 
method of having interaction terms appear in the 
energy and the three generators of the Lorentz 
group which give the three infinitesimal displace- 
ments in velocity. 


DYNAMICAL MODEL OF ELEMENTARY 
PARTICLES BASED ON AN ANALOGY WITH 
SUPERCONDUCTIVITY. I. Y. Nambu and 

G. Jona-Lasinio,* The Enrico Fermi Institute 
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for Nuclear Studies and the Department of Phys- 
ics, The University of Chicago, Chicago, Illinois 
(Received October 27, 1960). 


It is suggested that the nucleon mass arises 
largely as a self-energy of some primary ferm- 
ion field through the same mechanism as the 
appearance of an energy gap in the theory of 
superconductivity. The idea can be put into a 
mathematical formulation utilizing a generalized 
Hartree-Fock approximation which regards real 
nucleons as quasi-particle excitations. We con- 
sider a simplified model of nonlinear four - 
fermion interactions which allows a y,-gauge 
group. An interesting consequence of the sym- 
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metry is that there arise automatically pseudo- 
scalar zero-mass bound states of a nucleon-anti- 
nucleon pair which may be regarded as an 
idealized pion. In addition, massive bound 
states of nucleon number Zero and two are pre- 
dicted in a simple approximation. 

The theory contains two parameters which 
can be explicitly related to the observed nucleon 
mass and the pion-nucleon coupling constant. 
Some paradoxical aspects of the theory in con- 
nection with the y, transformation are discussed 
in detail. 

*Fulbright Fellow, on leave of absence from Istituto 


di Fisica dell’Universita, Roma, and Istituto Nazionale 
di Fisica Nucleare, Sezione di Roma, Italy. 








